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PREFACE 

This book was undertaken at the suggestion of my lamented 
colleague Professor Benjamin Osgood Peirce, and with the promise 
of his collaboration. His untimely death deprived me of his invalu- 
able assistance while the second chapter of the work was still 
unfinished, and I have been obliged to complete my task without 
the aid of his remarkably wide and accurate knowledge of Mathe- 
matical Physics. 

The books to which I am most indebted in preparing this treatise 
are Thomson and Tait's ^^ Treatise on Natural Philosophy," Watson 
and Burbury's "Generalized Coordinates," Clerk Maxwell's "Elec- 
tricity and Magnetism," E. J. Eouth's "Dynamics of a Rigid 
Body," A. G. Webster's " Dynamics," and E. B. Wilson's " Advanced 
Calculus." 

For their kindness in reading and criticizing my manuscript I am 
indebted to my friends Professor Arthur Gordon Webster, Professor 
Percy Bridgman, and Professor Harvey Newton Davis. 

W. E. BYERLY 
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CHAPTER I 

INTKODUCTION 

1. Coordinates of a Point. The position of a moving particle 
naay be given at any time by giving its rectangular coordinates 
a;, y, z referred to a set of rectangular axes fixed in space. It 
may be given equally well by giving the values of any three 
specified functions of x^ y, and 2, if from the values in question 
the corresponding values of x^ y^ and z may be obtained uniquely. 
These functions may be used as coordinates of the point, and 
the values of x^ y, and z expressed explicitly in terms of them 
serve as formulas for transformation from the rectangular sys- 
tem to the new system. Familiar examples are 'polar coordi- 
nates in a plane, and cylindrical and %pherical coordinates in 
space, the formulas for transformation of coordinates being 

respectively 

x = r cos 0, ] x = r cos ^, 



x = r cos <^, 
y = r sin 0, 



V 



(3) 



(1) y = r sin 0, i- (2) y = r sin Q cos 0, 

2 = 2, J 2 = r sin Q sin 0. 

It is clear that the number of possible systems of coordinates 
is unlimited. It is also clear that if the point is unrestricted in 
its motion, three coordinates are required to determine it. If it 
is restricted to moving in a plane, since that plane may be taken 
as one of the rectangular coordinate planes, two coordinates are 
required. 

The number of independent coordinates required to fix the 

position of a particle moving under any given conditions is 

called the number of degrees of freedom of the particle, and is 

1 



2 INTRODUCTION [Art. 2 

equal to the number of independent conditions required to fix 
the point. 

Obviously these coordinates must be numerous enough to fix 
the position without ambiguity and not so numerous as to render 
it impossible to change any one at pleasure without changing 
any of the others and without violating the restrictions of the 
problem. 

2. Dynamics of a Particle. Free Motion. The differential 
equations for the motion of a particle under any forces when 
we use rectangular coordinates are known to be 

mx — X, ] • 

mi/ = F, - (1)* 

niz — Z. 

X, F, and Z^ the components of the actual forces on the particle 
resolved parallel to the fixed rectangular axes, or rather their 
equivalents mx^ my^ mz^ are called the effective forces on the parti- 
cle. They are of course a set of forces mechanically equivalent 
to the actual forces acting on the particle. 

The equations of motion of the particle in terms of any other 
system of coordinates are easily obtained. 

Let g'j, q^y jg, be the coordinates in question. The appropriate 
formulas for transformation of coordinates express x^ y, and z in 
terms of q^^ q^, and q^. 

For the component velocity x we have 

dx , dx , ^ dx , 
dq^ ^' dq^ dq^ 

and x^ y, and z are explicit functions of j^, j^, j^, q^^ ^g, and j^, 
linear and homogeneous in terms of j^, q^^ and q^. 

* For time derivatives we shall use the Newtonian fluxion notation, so that 

we shall write x for — -, x for — -• 

at dt^ 



Chap. I] FEEE MOTION OF A PAETICLE 3 

We may note in passing that it follows from this ffU5t that 
a?, y^ and ^ are homogeneous quadratic functions of q^^ q^^ 

and q^. p- p 

Obviously —- = —-; (2) 

J . d dx S^x , , d^x . , d^x . 
and smce = — - a, H o^ H (7., 

A dx ___d^x , d^x . d^x . 

^^ = ^. (3) 

dtdq^ dq^ ^ ^ 

Let us find now an expression for the work iqW done by the 
effective forces when the coordinate q^ is changed by an infini- 
tesimal amount iq^ without changing q^ or q^. If 8a:, 8y, and 
hz are the changes thus produced in a:, y, and z, obviously 

iqW = m [^xSx + ySy + zSz] 

if expressed in rectangular coordinates. We need, however, to 
express BgWin terms of our coordinates q^, q^, and q^. 



-T ..dxd/. dx\ . d dx 

JN ow X — = — [x — ) — X ; 

dq dt\ dql dtdq 



but by (2) and (3) 



dx _dx , d dx _dx 

dq^ " a j, dt dq^ ~ dq^ 



„ ..dx_d(.dx\ .dx _d d (3?\ d /3?\ 

Hence x-^-^x-J-x----[^^J--\^-j; 

and therefore 8,.^= [| g - 1|] «?.' (4) 

■where r=^[i» + ^ + 2*] 

and is the kinetie energy of the particle. 



1^ = Q. (5) 
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To get our differential equation we have only to write the 
second member of (4) equal to the work done by the actual 
forces when q^ is changed by Bq^. 

If we represent the work in question by Qj^q^, our equation is 

dt dq^ cq^ 

and of course we get such an equation for every coordinate. 

It must be noted that usually equation (5) will contain q^ 
and q^ and their time derivatives as well as q^^ and therefore can- 
not be solved without the aid of the other equations of the set. 

In any concrete problem, T must be expressed in terms of q^^ 
q^y q^y and their time derivatives before we can form the expres- 
sion for the work done by the effective forces. Q^Bq^t Q^%^ 
QjSq^y the work done by the actual forces, must be obtained 
from direct examination of the problem. 

3. (a) As an example let us get the equations in polar coor- 
dinates for motion in a plane. 

Here x=r cos <^, !/ = f sin (f). 

x'-{-f = v^ = r^-\-r'4>% 

and T = ^[r' + t^c^^]. 



dT 



dr 
dT 



r = mr, 



= mr^^. 



dr 

i^W = m[r — r<^^] hr = Rhr 
if R is the impressed force resolved along the radius vector. 

— = 7wr^9, 



d4> 



= 0. 



(to 
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if 4> is the impressed force resolved perpendicular to the radius 
vector. 

In more f amihar form 



m 



r dt\ dt) 



(J) In cylindrical coordinates where 

x = r cos <^, y = r sin <^, « = «, 



T = 


= 2[^ + ^<^' + ^]. 




dT 
cr 












— — = mz. 



or m 



ai 

h^W = 771 [r — r<^^] 8r = Rir^ 
8JV = m'zSz = ZSz ; 

■ r dt\ dt) ' 
dh _ 
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(c} In spherical coordinates where 

x = r cos 0j i/ = rsm0 cos <^, 2; = r sin ^ sin <^9 

r = 5[r^ + 7^^H-7^sin^^<^^], 

ar 

^ = Tnr [^ + sin" ^<^n, 
— - = mr^r, 

-— -. = mr^ sin ^ cos 0<i>\ 

00 

— ;- = mr^ sin^ 6d>. 
S^W = mlr-r (p" + sin* ^<^')] hr = ESr, 
SeW= m\j (fd^ - r" sin ^ cos O^A B0 = erB0, 

S^W= m^(f sin' 0^^ h<f> = a)r sin 0h4> ; ' 






f[l(''f)-'"-'-Kf)>«. 



m 
rsin 



4. Dynamics of a Particle, Constrained Motion, If the particle 
is constrained to move on some given surface, any two inde- 
pendent specified functions of its rectangular coordinates a;, y, z, 
may be taken as its coordinates q^ and q^^ provided that by the 
equation of the given surface in rectangular coordinates and 
the equations formed by writing q^ and q^ equal to their values 
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in terms of a?, y, and z the last-named coordinates may be 
uniquely obtained as explicit functions of q^ and q^. For when 
this is done, the reasoning of Art. 2 will hold good. 

If the particle is constrained to move in a given path, any 
specified function of x^ y, and z may be taken as its coordinate 
g'j, provided that by the two rectangular equations of its path 
and the equation formed by writing q^ equal to its value in 
terms of x^ y, and z the last-named coordinates may be uniquely 
obtained as explicit functions of q^ For when this is done, the 
reasoning of Art. 2 will hold good. 

5. (a) For example, let a particle of mass w, constrained to 
move on a smooth horizontal circle of radius a, be given an 
initial velocity F", and let it be resisted by the air with a force 
proportional to the square of its velocity. 

Here we have one degree of freedom. Let us take as our 
coordinate q^ the angle 6 which the particle has described about 
the center of its path in the time t. 

and we have —^^mc^d. 

Our dijBferential equation is 

mame = - ha'd^W, 

" 1c Jin 

which reduces to -\ — aff^ = 0, 

m 



or 



at m 



Separating the variables, 

-^-\ — adt = 0, 

Integrating, ~-:H — cit= C = 

m V 
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1 


ma-\-k Vat 


de 


mV 
mV 


dt 


ma 4- k Vat 



tea 

and the problem of the motion is completely solved. 

(J) If, however, we are interested in i?, the pressure of the 
constraining curve, we must proceed somewhat differently. We 
have only to replace the constraint by a force B directed toward 
the center of the path. There are now two degrees of freedom, 
and we shall take Q and the radius vector r as our coordinates 
and form two differential equations of motion. 





^-if 


^' + r'O'l 




dT 
d6 


mr^d, 




dT 
dr 


mr, 




dT 
dr 


mrO^. 


d 

m-z- 

dt 


(fd) 80 = 


- kr'd^Bd. 



m 
To these we may add 



(1) 
(r - rd^) Sr = - Rhr. (2) 



r = a. 



Whence d* -h — ^ = 0, as before, (3) 

m ^ '^ 

and B = ma^. (4) 
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(<?) Let us now suppose that the constraining circle is rough. 
Here, since the friction is /a (the coefficient of friction) multi- 
plied by the normal pressure, B will be needed, and we must 
replace the constraint by ^ as before. 

We have now 

at 





m\r — r^] hr = — Rir^ 


and 


r = a. 


Whence 


B = 7nad^^ as before. 


• 


m ma 




_m 



Replacing — in Art. 5, (a), (1), by — H- /i, 
m m 



we have <,_^,c«[l + (^ + .)0]. (1) 

m 

EXAHPLES 

1. Obtain the familiar equation — ^-|--sin^ = for the 

1 J 1 at a 

sunple pendulum. 

2. Find the tension of the string in the simple pendulum. 

An%. B=^m\gQ.Q%Q-\-a\—\ . 

3. Obtain the equations of the spherical pendulum in terms 
of the spherical coordinates 6 and <^. 

An%. -Bine cos0<}>^-^^8m = 0. sin^^<^=a 

a 

6. (a) The constraint may not be so simple as that imposed 
by compelling the moving particle to remain on a given surface 
or on a given curve. 
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Take, for example, the tractrix problem, when the particle 
moves on a smooth horizontal plane. 

Let a particle of mass m, attached to a string of length a, rest 
on a smooth horizontal platie. The string lies straight on the 
plane at the start, and then the end not attached to the particle 
is drawn with uniform velocity along a straight line perpendicu- 
lar to the initial position of the string and lying in the plane. 
Let us take as our coordinates 2%, the distance traveled by 
that end of the string which is not attached to the particle, 
and 6^ the angle made by the string with its initial position. 
Let B be the tension of the string and n the velocity with 
which the end of the string is drawn along. Let X, F, be the 
rectangular coordinates of the particle, referred to the fixed 
line and to the initial position of the 
string as axes. 

X=x — a sin ^, 

F = a cos 6 ; 

X= X — a cos dd^ 

r=-asin^^. 




O 



X 



r=^(X^-hr^)=^[2:«H-a^^-2acos^i^]. 



dT 



ex 



dT 



r = m[^x — a cos ^^], 



--. = m [a^^ — a cos ffx]j 



dd 



= ma sm 



in 6xd, 



d^. 



m 



m — \x — a cos ^^] Sa: = -R sin Ohx. 
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Adding the condition x = nt^ 

and reducing, — ma [cos 06 — sin 6^'\ = ^ sin ft 

ma^'e = 0. 

R = maff^. 
Integrating, ^ = (7 = -. 

a 

The particle revolves with uniform angular velocity about 
the moving center, and the pull on the string is constant. 

(J) A particle is at rest in a smooth horizontal tube. The 
tube is then made to revolve in a horizontal plane with uniform 
angular velocity o). Find the motion of the particle. 

Suggestion. Take the polar coordinates r, <^, of the particle as 
our coordinates, and let R be the pressure of the particle on 
the tube. 

dT 

-^ = rnr, 

or 

^^ ^JL 
— ^ = rmr<p. 

m [r — r(f>^'\ Br — 0. 
m ■— (r^<^) 8<^ = ErS<f>, 

fA/tr 

Adding the condition <^ = ©^ 

and reducing, r — ©V = 0, 

2 mcorr = Er. 
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Solving, r = -4 cosh (Dt'\-B sinh ©^ 

r = a and f = at the start. 

Hence r = a cosh a)t = a cosh (f>, 

E = 2 mcuD^ sinh ©^ = 2 wa©^ sinh <^. 

If we are interested only in the motions and not in the reac- 
tions, problems (a) and (J) can be solved more simply. If in 
each we were to use one less coordinate, only in (a) and r 
only in (J), rectangular coordinates JT, F, for the particle could 
be obtained whenever the time was given, and therefore could 
be expressed explicitly in terms of ^ or r and t, A careful 
examination of Art. 2 will show that the reasoning is extended 
easily to such a case, and that the work done by the effective 

forces when g^ only is changed is still — t-, — -r— ^Q.» It i 

to be noted, however, that when the rectangular coordinates are 
functions of ^ as well as of q^, q^, etc., the energy T is no longer 
a homogeneous quadratic in q^, q^j etc. 

(a') jr=w^ — asin^, 

F=acos^; 

r=-asin^A 



IS 



—-.=^m \a^d — an cos 6\ 
dd 

-— r = rtian sm &a. 
dd 

m -j(«^^ — «^ cos ^) — an sin ^^ S^ = 0, 

^ = - , as before. 
a 
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(6') r = |[r^ + a,V]. 

dT 
Cr 

dT 2 

— — = nmr, 

dr 

w [r — 6)V] Sr = 0. 
r = a cosh o)^, as before. 

EXAMPLES 

1. A particle rests on a smooth horizontal whirling table and 
is attached by a string of length a to a point fixed in the table 
at a distance b from the center. The particle, the point, and 
the center are initially in the same straight line. The table is 
then made to rotate with uniform angular velocity ©• Find the 
motion of the particle. 

Suggestion, Take as the single coordinate the angle made 
by the string with the radius of the point. Let -ST, Y, be the 
rectangular coordinates of the particle, referred to the line ini- 
tially joining it with the center and to a perpendicular thereto 
through the center as axes. 

Then X=h cos a)t'\-a cos (^ -|- (of)^ 

and F= J sin ©^ -}- a sin (^ -}- (of). 

r=5[JV + a2(a) + ^y-h2a6Q)COS^(a)-|-^)], 

and H sin ^ = ; 

a 

and the relative motion on the table is simple pendulum motion, 

the length of the equivalent pendulum being I = 7-^ • 

2. A particle is attracted toward a fixed point in a horizontal 
whirling table with a force proportional to the distance. It is 
initially at rest at the center. The table is then made to rotate 
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with uniform angular velocity od. Find the path traced on the 
table by the particle. 

Suggeatioru Take as coordinates a;, y, rectangular coordinates 
referred to the moving radius of the fixed point as axis of 
abscissas and to the center of the table as origin. Let X, F, be 
the rectangular coordinates referred to fixed axes coinciding 
with the initial positions of the moving axes. 

X= X cos (ot — y sin <ot^ F= x sin ©^ -f y cos w^. 

Whence come tw [^ — 2 o)^ — w^a;] = — /i (a; — a), 

m\_y '\- 2 cax — co^y^ = — fiy. 

If ©^ = — » the solution is easy and interesting. 
m 

x — 2a)y = aG)% (1) 

y-{-2a)x=0. (2) 

Integrating (2), y -\- 2 (ox = 0, 

Substituting in (1), a; -j- 4 oy^x = aG)\ (3) 

Multiplying (3) by 2 a?, and integrating, 

d^-\-4: © V = 2 aoy^x. 



Whence 



x= 2a)-v 


2 ^- 


o J. -i4a; 
2a)t = vers ^ — . 

a 


a? = J [1 — cos 2 G)t'], 


a j.^ 


lG)t — sm2 cot]. 



Replacing 2 ©^ by ^, a; = j [1 — cos ^], 
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and the curve traced on the table is the cycloid generated by a 
circle of radius j rolling backward along the moving axis of Y. 

7. A System of Particles. If instead of a single particle we 
have a system of particles, free, or connected or otherwise con- 
strained, Twai, my^ mZy are the effective forces on the particle P. 
The effective forces on all the particles are spoken of as the 
effective forces on the system and are mechanically equivalent 
to the set of actual forces on the system. 

T, the kinetic energy of the system, is the sum of the kinetic 
energies of all the particles. 

T=^^[f + f + i?^. 

If hW is the work done by the effective forces in any supposed 
infinitesimal displacement of the system, 

hW= 2w \xix -I- yhy -\- zhz\. 

If the particles of a moving system are subjected to connec- 
tions or constraints, these connections or constraints may or 
may not vary with the time. In the latter case a set of any n 
independent variables q^^ q^, • • •, q^j such that when they and 
the connections and constraints are given, the position of every 
particle of the system is uniquely determined, and such that 
when the positions of all the particles of the system are given, 
9v ?2» • • •' ?»' follow uniquely, may be taken as coordinates of 
the system ; and n is called the number of degrees of freedom 
of the system. 

In the former case a set of variables g^^, j^, • • •, g,, such that 
when they and the time are given, the position of every particle 
of the system is uniquely determined, and such that when the 
positions of all the particles of the system and the time are 
given, q^^ y^, . . ., j^, follow uniquely, may be taken as the coor- 
dinates of the system ; and n is called the number of degrees of 
freedom of the system. 
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The equations expressing the connections and constraints in 
terms of the rectangular coordinates of the particles and of the 
coordinates y^, q^, • • •, j^, of the system are often called the 
geometrical eqiuitiona of the system and may or may not contain 
the time explicitly. In the latter case the geometrical equations 
make it possible to express the coordinates x, y, 2, of every 
point of the system explicitly as functions of the g-'s; in the 
former case, as functions of t and the j's. 

The geometrical equations must not contain explicitly either 
the time derivatives of the rectangular coordinates of the parti- 
cles or those of the coordinates q^, q^, • • •, q^, of the system 
unless they can be freed from these derivatives by integration. 

Examples of geometrical equations not containing the time 
explicitly are the formulas for transformation of coordinates in 
Arts. 1 and 3, and the equations for X and Y in Art. 6, (a). 

Geometrical equations containing the time are the equations 
for X and Y in Art. 6, (a'), and in Art. 6, Exs. 1 and 2. 

The work, S^ PT, done by the effective forces when q^ is 
changed by Bq^ without changing the other g-'s is proved to be 

by reasoning similar to that used in Art. 2. For the sake of 
variety we take the case where the geometrical equations 
involve the time. 

Here x =f\t, q^, q^, . . ., gj. 



dx dx . dx , ^ dx , 

dt dq^^^ ag/2 dq^^ 

and is an explicit function of ^, q^, q^, • . •, q^, q^, q^, • • ., q^. 

dx ^ dx 

, . d dx d^x , d\ , , d^x , , ^ d^x , 

and smce 1: ^ = r;7-H-7-i?i + ^-r"92 "• '■r"F~?»' 



(1) 
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f«» 



X 



dx ^x , d^x . , ^x . , , 8'x 
and — ^ 1 g H g + • • • H 

— [—\ = — r2'J 

dt\dqj dq; ^^ 

K A HJ dt\dqj dt\ dqj dq^^"-^ ^^ 

'^ dq- dt[dq\2)\ dqX^y 

and therefore h,W= [|(g) - 1|] 89,, 

and if Q^q^ is the work done by the actual forces when q^ is 
changed by hq^, d dT dT 

If the geometrical equations do not contain the time, the same 
result is seen to hold, and in this case it is. to be noted that 
since x is homogeneous of the first degree in the time derivatives 
of the coordinates, that is, in ^^, q^^ . • ., 5-,, the kinetic energy T 
is a homogeneous quadratic in q^^ q^, • • ., q^. 

Generally every one of the n equations of which equation (3) 
is the type will contain all the n coordinates 5'^, q^^ • • •, q^, and 
their time derivatives, and can be solved only by aid of the 
others. That is, we shall have n coordinates and the time con- 
nected by n simultaneous differential equations no one of which 
can ordinarily be solved by itself. 

If the forces exerted by the connections and constraints do 
no work, they will not appear in our differential equations. 
Should we care to investigate any of them, we have only to 
suppose the Constraints in question removed and the number 
of degrees of freedom correspondingly increased, and then to 
replace the constraints by the forces they exert and to form 
the full set of equations on the new hypothesis. 
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8. A System of Particles. Illustrative Examples, (a) A rough 
plank 10 feet long rests pointing downward on a smooth plane 
inclined at an angle of 30° to the horizon. A dog weighing 
aH much as the plank runs down the plank just fast enough 
to ke(ip it from slipping. What is his velocity when he reaches 
itH lower end? 

Hero we have two degrees of freedom. Take x^ the distance 
of tli(} upper end of the plank from a fixed horizontal line in 
the plane, and y, the distance of the dog from the upper end 
of tlui plank, as coordinates, and let R be the backward force 
exerted by the dog on the plank, and m the weight of the dog. 

dT 

— = «t[2i + y], 

dT 

— = m[x-\-y-]. 

w [2 ij + y] 8a; = 2 w^ sin 30° ix, 
7n [i + jf] % = [-B -f mg sin 30°] hy. 

2x + y=:g, 

,. .. R g 
ml 

a: = a constant, 

f^^'^gy^c^^gy, 

y=.y/2gy. 
y = 32, nearly. 

R , q 
ml 



Hy hypothesis, 
and therefore 



When y = 16, 
Since 



i^ = 



mg 
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(6) A weight 4w is attached to a string which passes over 
a smooth fixed pulley. The other end of the string is fastened 
to a smooth pulley of weight m, over which passes a second 
string attached to weights m and 2 m. 
The system starts from rest. Find the 
motion of the weight 47w. 

Two coordinates, x^ the distance of 4 m 
below the center of the fixed pulley, and 
y, the distance oi 2 m below the center 
of the movable pulley, will suffice. The 
velocities are 



Am 




r^ 




m 



1 2m 



] 



X for 4 m, 
— X for movable pulley, 

— a; + y for 2 w, 

— x — y for 7w. 

T= ^[imd? -hmx" + 2m(y -^ xy -{-mOb-^- yy] 

dT 

— = m[8i-y], 

dT 

— = m[^y-x]. 

?w [8 ai — y] Sa; = [4 ?w^ — mg — 2 mg — mg"] Sx, 

7W [3 y — i] Sy = [2 mg — mg"] Sy, 

8x-y = 0. 

y-x = g. 

2Sx = g. 
q 
23 

The weight 4 m will descend with uniform acceleration equal 
to one twenty-third the acceleration of gravity. 

((?) The dumh-hell problem. Two equal particles, each of mass 
7n, connected by a weightless rigid bar of length a, are set 
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• 

moving in any way on a smooth horizontal plane. Find the 
»ubH(5(][uent motion. 

Wo have tliree degrees of freedom. Let a:, y, be the rectangu- 
lar coordinates of the middle of the bar^ and the angle made 
by tlie bar with the axis of X. 

The rectangular coordinates of the two particles are 



(a;— -coH^, y — -mi0\ and (a;-f--cos^, y-f--sintfj; 



tlioir vclociticH are 



and -Mi - ^ sin edj+ U + ^ cos d^Jl 

7'= ^1 ir" + / + ^ #' +.(« sin ^ - a co8d)^(i: + y) 

+ z' + f + ^ ^ + (^a coad - a sme)&(x + ^)\ 

dT - . 
—- = 2mx, 

ox 

dT _ . 
dT ma' x 

2 mxSx = 0, 
2 mi/Si/ = 0, 

^680 = 0. 
a; = 0, 

y = o, 

6' = 0. 
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Hence the middle of the bar describes a straight line with 
uniform velocity, and the bar rotates with uniform angular 
velocity about its moving middle point. 

EXAMPLE 

Two Alpine climbers are roped together. One slips over a 
precipice, dragging the other after him. Find their motion 
while falling. 

Ans. Their center of gravity describes a parabola. The rope 
rotates with uniform angular velocity about their moving center 
of gravity. 

Qd) Two equal particles are connected by a string which 
passes through a hole in a smooth horizontal table. The first 
particle is set moving on the table, at right angles with the 
string, with velocity Vo^ where a is the distance of the particle 
from the hole. The hanging particle is drawn a short distance 
downward and then released. Find approximately the subse- 
quent motion of the suspended particle. 

Let X be the distance of the second particle below its position 
of equilibrium at the time f, and the angle described about 
the hole in the time t by the first particle. 

T^'^ld^ + d^ + Ca^xyd^. 

ox 

--— = — 7w (a — x)0*j 

ex 

— - = m(a — xyd, 
cd 

m\^x-\-(a-'X) 6^'] 8x = mgSaCj (1) 

m^Ka-xyd'jBe^O. (2) 

2x+(a-x)^ = g. (3) 

(a - xf6 = C = a V^, (4) 
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since (2) holds good while the hanging particle is being drawn 
down as well as after it has been released. 






2x-\ ^ = 0, approximately, 

a 

and i -h ^ a; = 0. 

2a 

For small oscillations of a simple pendulum of length ?, 

+ ^0 = 0. 

Therefore the suspended particle will execute small oscilla- 
tions, the length of the equivalent simple pendulum being ^ a. 

EXAMPLE 

A golf ball weighing one ounce and attached to a strong 
string is " teed up " on a large, smooth, horizontal table. The 
string is passed through a hole in the table, 10 feet from the 
ball, and fastened to a hundred-pound weight which rests on a 
prop just below the hole. The ball is then driven horizontally, 
at right angles with the string, with an initial velocity of a 
hundred feet a second, and the prop on which the weight rests 
is knocked away. 

(a) How high must the table be to prevent the weight from 
f alUng to the ground ? (6) What is the greatest velocity the 
golf ball will acquire ? Ans. (a) 8.96 ft. (V) 963.4 ft. per sec. 

9. Rigid Bodies. Two-dimensional Motion. If the particles of 
a system are so connected that they form a rigid body or a 
system of rigid bodies, the reasoning and formulas of Art. 7 
still hold good. 
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(a) Let any rigid body containing a horizontal axis fixed in 
the body and fixed in space move under gravity. Suppose that 
the body cannot slide along the axis. Then the motion is 
obviously rotational, and there is but one degree of freedom. 
Take as the single coordinate the angle made by a plane con- 
taining the axis and the center of gravity of the body with a 
vertical plane through the axis. 

Let h be the distance of the center of gravity from the axis, 
and k the radius of gyration of the body about a horizontal 
axis through the center of gravity. Then 

T = |(A^ + A2)^^ (v.App.A,§§5andlO) 

m (A' + le) 080 = - mgh sin 0S0. 

and we have simple pendulum motion, the length of the equiv- 
alent Simple pendulum being 

h' + It" 



1 = 



h 



(^) Two equal straight rods are connected by two equal 
strings of length a fastened to the ends of the rods, the whole 
forming a quadrilateral which is then suspended from a hori- 
zontal axis through the middle of the upper rod. The system 
is set 'moving in a vertical plane. Find the motion. 

Take as coordinates <f>, the inclination of the upper rod to 
the horizon, and 0, the angle made with the vertical by a line 
joining the point of suspension with the middle of the lower 
rod. From the nature of the connection the rods are always 
parallel 

Let k be the radius of gyration of each rod about its center 
of gravity. 
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(V. App. A, § 10) 

— - = mcrQ. 

2 mlc'^h^ = 0, 
^ = 0. 

a 

and the rods revolve with uniform angular velocity while the 
middle point of the lower rod is oscillating as if it were the 
bob of a simple pendulum of length a. 

(f) If an inclined plane is just rough enough to insure the 
rolling of a homogeneous cylinder, show that a thin hollow 
drum will roll and slip, the rate of slipping at any instant being 
one half the linear velocity. 

Let X be the distance the axis of the cylinder has moved 
down the incline, the angle through which the cylinder has 
rotated, a the radius of the cylinder, and a the inclination of 
the plane. Call the force of friction F. 

r = 5 [i^ + A^^]. 

dT 

ox 

mxix = [mg sin a — F'] hx^ 
mi^eSd = FaSd. 
If there is no slipping, x = aOy 

rnx =^ mg sma " F^ 
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■ • 

ml^ - = F(u 
a 

Hence -rr^ = f^g sin a—F^ 

For a solid cylinder, A^ = — , 

F= ^ Tw^sinor; 
R = mg cos a, 
where R is the pressure on the plane; 

F 1^ 
^ = - = - tan a, 

where fi is the coefiBcient of friction. 
For a hollow drum, A = a, 

i^' = 1 7w^ sin a, 

— = - tan a, 
^ 2 

ft < ^ tan a, 
and the drum will slip. 

For the drum, then, 

F = fiE = fimg cos (x = ^ mg sin a, 
mx = mg sin a: — ^ mg sin a = |^ mg sin a, 
x = ^gt sin or. 

mc?0 = Fa = -mgmi a, 

o 

a^ = J (7 sin a, 

S=x — au=^- mgt sin a = - , 
where 8 is the rate of slipping. 
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EXAMPLES 

1. A sphere rotating about a horizontal axis is placed on a 
perfectly rough horizontal plane and rolls along in a straight 
line. Show that after the start friction exerts no force. 

2. A sphere starting from rest moves down a rough inclined 
plane. Find the motion, (a) What must the coefficient of 
friction be to prevent slipping ? (6) If there is slipping, what 
is its velocity? 

An%. (a) ft > I tan a. (6) S = gt[^ma — ^ficma'\. 

3. A wedge of mass Jlf having a smooth face and a perfectly 
rough face, making with each other an angle a, is placed with its 
smooth face on a horizontal table, and a sphere of mass m and 
radius a is placed on the wedge and rolls down. Find the motion. 

Let X be the distance the wedge moves on the table, and y 
the distance the sphere rolls down the plane. 

Note that T==^\^M-^m']a?-^'^Y-^ 

Arts. (m + M^x — mt/ cos a = 0, ^y — xcosa = ^gf sma. 

10. Rigid Bodies. Three-dimensional Motion, (a) A homo- 
geneous sphere is set rolling in any way on a perfectly rough 
horizontal plane. Find the subsequent motion. 

Let X, y, a, be the coordinates of the center of the sphere 
referred to a set of rectangular axes fixed in space; two of 
which, the axes of X and F, lie in the given plane. Let OAy 
OB, OCy be rectangular axes fixed in the sphere and passing 
through its center ; let OX, Y, OZ, be rectangular axes through 
the center of the sphere parallel to the axes fixed in space; and 
let 6, ^, y^ be the Euler's angles (v. App. A, § 8). Take x, y, 
6, <f>, and yjr as our coordinates. The only force we have to 
consider is F, the friction, and we shall let F^ and F^ be its 
components parallel to the axes OX, OY, respectively. 



T=:^ 



|[i« + ^ + FK + a,J + a,^)], 
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where ©^ = — ^ sin -^/r + ^ sin eos -^j 

ft)y = ^ eos '^ + <j> sin ^ sin y^, 

<o^ = <j> cos -\- yjr. (v. App. A, § 8) 

Hence T=^[i;^ + y' + **(^ + <^* + ^*+2co8^<^)]. 

We get »ni = F^, (1) 

% = F^, (2) 

nt*» ^ [^^ + cos ^<^] = 0, (3) 

ft 

mT^ --- [<^ + cos 0'>^'] = — aF^ sin ^ sin -^/r 4- a^^ sin cos -^^ (4) 

• • • • 

mlr [0 + sin 0^y^'\ = — aF^ cos '^ — ai'y sin -^ ; (5) 

and as there is no slipping, 

X — aft)y = x — a(6 cos -^ + <^ sin ^ sin -^/r) = 0, (6) 

y + a®^ = y + a (— (? sin -^ + ^ sin ^ cos '^) = 0. (7) 
From (4) and (5), 

mJ^ [sin ^r— QJ>-\- cos ^'<^) + sin ^ cos -^ (^ + sin 0<f>yjr^'\ 

= -aF,sm0, (8) 

flu • • •• • • 

wAr^ [cos '^ T (<A + c^s ^'^) — sin tf sin -^/r (^ -f- sin ^<^) ] 

= aF^ sin ^. (9) 

Expanding the first members of (8) and (9) and eliminat- 
ing ^ by the aid of (3), we get 

mJ(^ [cos yjr0 — sin slrd-^ + sin sin yjrff) + cos sin -^/r^c^ 

4- sin ^ cos yjr^'] = — ai^^, (10) 

mlc^ [— sin '^/rd' — cos sjrd^ 4- sin ^ cos yfr^ + cos cos yjr^^ 

— sin ^ sin yfr^'l = ai^^,. (11) 
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But the first members of (10) and (11) are obviously mlc^ —-^ 

J dt 

,»d(o_. .. , ^^ , .^. , ... 



dt 


ny» xaeuue, uy ^u^ 




x = -'aF^, 

a 






Substituting in (1) 


and (2), we get 




a: = — max^ 

a 




y= may\ 


whence 


a;=0, 




y = 0. 


From (7) and (6), 


^-^' 




^=^^ 


and from (3), 


dt 


Finally, 


F, = 0, 




F.. = 0. 



V 



Hence the center of the sphere moves in a straight line with 
uniform velocity, and the sphere rotates with uniform angular 
velocity about an instantaneous axis which does not change its 
direction, and no friction is brought into play after the rolling 
begins. 

(6) The billiard ball. Suppose the horizontal table in (a) is 
imperfectly rough, coefficient of friction /i, and suppose the 
ball to slip. 

Take the same coordinates as before, and equations (1), (2), 
(3), (4), (5), (8), (9), (10), and (11) still hold good. Let a 
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be the angle the direction of the resultant friction, F=fimg^ 
makes with the axis of X, and let S be the velocity of slip- 
ping, that is, the velocity with which the lowest point of the 
ball moves along the table. Of course the directions of F and 8 
are opposite. 

Let S^ and S^ be the components of S parallel to the axes of 
X and F. We have 

— S cos a= S^ = x — aft)y, 
and — aS^ sin a = /Sy = y -f- <3ta)^. 

F^ = fimg cos a, 
and Fy = fimg sin a. 

dS^ „ . da dS .. dto^ 

-rr- = iS sm Of -- — cos a --- = a; — a —-»• 

at dt dt dt 

dS„ ^ da . dS .. ^ dto^ 

-— ^ = — /Scosa-- — sma-— - = y + a-— 2. 

dt dt dt ^ dt 

From (1), x = fJLg cos a, 

and from (10), — * = — — /i^ cos a. 

XT „ , da dS a^ + Jc^ .-on 

Hence S^ma-- — cosa— = — — — ft^cosa, (12) 

CLZ CbZ fC 

and from (2) and (11), 

^ da . dS a^-j-Ii^ . .^on 

— Scoaa -- — sma — = — z^ — figsma. (13) 

CtZ CLZ rC 

Multiplying (12) by sin a and (13) by cos a, and subtracting, 

^1 = 0. (14) 

Multiplying (12) by cos a and (13) by sin a, and adding, 

dS a^ + Zc" 



dt y 



ng. (15) 
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Integrating (15), 5=^^-?l±^^ (16) \ 

From (14), « = ^o' 

and the direction of slipping does not change. 

If the axes are so chosen that the axis of X has the direction 
opposite to the direction of slipping, a = 0. Then 

y = o. 

These equations are familiar in the theory of projectiles, 
and the path traced on the table is a parabola so long as 
slipping lasts. 

Should y^ happen to be zero, the path degenerates into a 
straight line. 

When slipping stops, 

X — aft)y = 0, and y + a©^ = 0, 
aad we have the case treated in (a). 

EXAMPLE 

A homogeneous sphere is set rolling on a perfectly rough 
inclined plane. Find the path traced on the plane. 

Ans, A parabola. 

((?) 77ie ffyroscope. Suppose a rigid body containing a fixed 
point and having two of its moments of inertia about its prin- 
cipal axes through the fixed point equal. Obtain the differ- 
ential equations for its motion under gravity. 

We shall use Euler's angles with a vertical axis of Z. 

We have ©^ = ^ sin <^ — -^^ sin ^ cos <^, 

(o^ = 6 cos <^ -f- '^ sin ^ sin <^, 
©g = '^ cos 6 -{-f^. (v. App. A, § 8) 

r= I lAa>l+Aa>l + C^o^-] (v. App. A, § 10) 

= K-^ (J^ + s^^' ^^') + ^(^ cos ^ + <^)']. 
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^- = CQ<}r COB + 4>}, 
d<l> 

^'=AsiD!'^^jr-\'C COS (^jr COS + <^), 

d0 

-;~=A%m0 COS 0'^^ — Csia0(yjr cos -f- ^)'^. * 
du 

Our equations are C— (-^ cos ^ + ^) = 0, (1) 

|-[^sin''d^ + Cco8^(-^co8d + ^)] = 0, (2) 

AS — Asia cos 0-^^ +C am 0(^jr cos + ^^-^== mga sin tf. (3) 

From (1), yjrcoa0 + ^ = a, (4) 

where a: is the initial velocity about the axis of unequal moment. 

A sin' 0-^ + Ca cos 0=:L. (5) 

AS — Aain cos ^i^' + Ca sin ^-^^ = mga sin ^, (6) 

or substituting yjr from (5), 

v^ (X — Cacos^)(Xcos^--Ca) , . ^ ... 

-^4^ = ^^ r^Vs + ^« sm ^. (7) 

(tf) Obtain Euler's equations for a rigid body containing a 
fixed point. 

Here T = i [Acj^ + ^o),' + (7©«]. (v. App. A, § 10) 

O/TT 

— =^G)j [^ cos ^ + -^ sin ^ sin </>] 
^ +-^®2 [— ^ sin <^ + '^ sin ^ cos <^] 

= A(o^(o^ — B(ojX)^. 
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Whence [^^ - (^ ^ ^) ^i®,] ^ = ^^^ 
where N is the moment of the unpressed forces about the C axis. 

The remaining two Eulers equations follow at once from 
this by considerations of symmetry. 

11. In Arts. 2 and 7 it was shown that under slight limi- 
tations the coordinates of a moving particle or of a moving sys- 
tem coujd be taken practically at pleasure, and the differential 
equations of motion could be obtained by the application of a 
single formula. It does not follow, however, that when it comes 
to solving a concrete problem completely, the choice of coordi- 
nates is a matter of indifference. Different possible choices 
may lead to differential equations differing greatly in compli- 
cation, and as a matter of fact in the illustrative problems of 
the present chapter the coordinates have been selected with care 
and judgment. That this care, while convenient, is not essential 
may be worth illustrating by a practical example, and we shall 
consider the simple familiar case of a projectile in vacuo. 

Altogether the simplest coordinates are x and t/j rectangular 
coordinates referred to a horizontal axis of X and a vertical 
axis of Y through the point of projection. 



We have 


^=f(^ + 2^), 




dT 

ex 




dT 




mx= 0, 




mi/ = — mg. 
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Solving, X = v^, 

X = vj,. 

Let us now try a perfectly crazy set of coordinates, q. and q, 

where . , , 

q^ = x + iaxr^y, 

and J =a; — tan~*y. 

Proceeding in our regular way, we have 

* 2 

y = 2 sec' 2i^ ( j, - j,). 

^ = f [?? + 2 ?A + ?H sec^ 2lZ^^ (^J - 2 q^q^ + y^)] . 

g^ = 4 [?i + ?2 + sec* ^J-2^ (?i - ?,) J, 
^ _ ^ sec* 2lJl9? tan 2l:^^ ro _oV 

g^ = 4 [?i + 92 - sec" ^i-2^ (?, - ?,) J, 



^9 

^^ _ m ^^r-i gi - 92 

dq~ 4 2 "" 2 



= - ^ sec- ^--:.-- tan 2»-_^ ( j^ - ^^)«. 



f[9x + 92 + sec"2i^(^^-^^) + 8ec"2i_&tan2L_&G,-j,y] 



= _^sec«2i^. (1) 
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= ^8ec«2^. (2) 

Adding (1) and (2), |(9, + 9,) = 0. 

Whence ?i + ?, = ^ ^a 

and ?i + ?, = 2 v/. (3) 



VY nence 9i'rq^ = ^ r^, 

id ?i + ?a = 2 v/. 

Subtracting (2) from (1), 

ifsec' 2lZ^ (?,-?.) + sec* ^ tan 2i^ G, - j,)'] 

= — wi^ sec' 2i-— 2i . 

4 

Multiplying by — (j — y ), and integrating, 

m 

sec* 2i^=^ ( j, - j,)* = - 8 tan 2lZL^ + 4 ^. 

Bec« 2^ (y, - 9,) = 2 >Jt^, - 2^ tan 2i^ 

Let ^ = ^- 

sec' 2 -r = ^''» — 2 ^ tan «, 

,, sec'zrfz 

at = — » 

vv^ — 2^tanz 
« = -[v,-^y,'-2<7tanz] 



= ^[''»-\|<-2^tan2i£^ 
^^ 2 •'2 



V-'^-f'.'-fl 
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But from (3), il±^ = vj. 

Hence q^ = vjt + tan"* V "" 9 ' 

9, = t;,<-tan-*[v-^l- 

Of course this should agree with our first answer, and a 
moment's consideration shows that it does. 



We have Si±3i ^ ^c =- vj, 



gt* 



tan ^* "^ = y :=! vjt — —- y as before. 

12. The parameters j^, ^j, • • • that we have been using to 
fix the position of our moving particle or moving system are 
called generalized coordinates. Following the analogy of rectan- 
gular coordinates, the time derivative q^ of any generalized 
coordinate q^ is called the generalized component of velocity 
corresponding to qj^. It may be a linear velocity, or an angular 
velocity as in many of our problems, or it may be much more 
complicated than either as in our latest example. 

The kinetic energy T expressed in terms of the generalized 
coordinates and the generalized velocities is called the Lagran- 
gian expression for the kinetic energy. 

If we are using rectangular coordinates and dealing with a 
moving particle, 

dT 

ox 

and is the x component of the momentum of the particle. 

dT 
Following this analogy, — is called the generalized com- 

ponent of the momentum of the system, corresponding to the 
coordinate j^.. It is frequently represented by p^y and may be a 
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momentum, or a moment of momentum as in many of our prob- 
lems, or it may be much more complicated than either as in 
our latest example. 
Equations of the type 

are practically what are called the Lcigrangian equations of 
motion, although strictly speaking the regulation form of the 
Lagrangian equations is a little more compact and will be 
given later, in Chapter IV. 

Q^^ defined through the property that Qj^iqjg is the work done 
by the actual forces when 9^ is changed by Sj^, is called the 
generalized component oi force corresponding to 9^. It may be 
a force, or the moment of a force as in many of our problems, 
or it may be much more complicated than either as in our 
latest example. 

13. Summary of Chapter I. If a moving system has a finite 
number n of degrees of freedom (v. Art 7) and n independent 
generalized coordinates j^, q^^ • • •, y„, are chosen, the kinetic 
energy T can be expressed in terms of the coordinates and 
the generalized velocities 9^, y^, • • •, q^^ and when so expressed 
will be a quadratic in the velocities, a homogeneous quad- 
ratic if the geometrical equatio7i% (v. Art. 7) do not contain the 
time explicitly. 

The work done by the effective forces in a hjrpothetical infin- 
itesimal displacement of the system due to an infinitesimal 
change dqj^ in a single coordinate q^ is 



[dt dq, ~dq,\ ?*• 



If this is written equal to ©^^Sj^, the work done by the actual 
forces in the displacement in question, there will result the 
Lagrangian equation , ^^ ^^ 

dt dqj^ dqj^ 



- — =0*. 
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The n equations of which this is the type form a set of simul- 
taneous differential equations of the second order, connecting 
the n generalized coordinates with the time. When the complete 
solution of this set of equations has been obtained, the problem 
of the motion of the system is solved completely. 

It must be kept in mind that in order to obtain the value of 
a single coordinate or of a set of coordinates less in number 
than n it is generally necessary to form and to solve the com- 
plete set of n differential equations. 

We shall see, however, in the next chapter, that in certain 
important classes of problems some of these equations need not 
be formed, and that some of the coordinates can be safely 
ignored without interfering with our obtaining the values of the 
remaining coordinates; that, indeed, we may be able to handle 
satisfactorily some problems concerning moving systems having 
an infinite number of degrees of freedom. 



CHAPTER II 

THE HAMILTONIAN EQUATIONS. ROUTH'S MODIEIED 
LAGRANGIAN EXPRESSION. IGNORATION OF 

COORPINATES 

14. The Hamiltonian Equations. If the geometrical equations 
of the system (v. Art. 7) do not contain the time expUcitly 
and the kinetic energy T is therefore a homogeneous quad- 
ratic in 9j, 9j, • • •, ^,, the generalized component velocities, 
Lagrange s equations can be replaced by a set known as the 
Hamiltonian equations. 

The Lagrangian expression for the kinetic energy we shall 
now represent by T^. 

Let « =— -^, « =-— ^, etc. be the generalized component 

dq^ dq^ 

momenta. Then jt?^, jc^j, • • • are homogeneous of the first degree 
in 9j, 9a, • • •. Express q^, J^, • • • in terms of p^,p^, • • •» 9i» 92» • • •» 
noting that they are homogeneous of the first degree in terms 
of jt?j, /?j» • • •» and substitute these values for them in T^, which 
will thus become an explicit function of the momenta and the 
coordinates, homogeneous of the second degree in terms of the 
former. This function is called the Hamiltonian expression for 
the kinetic energy, and we shall represent it by T^. Of course 

Ti = Tp. (1) 

^ By Euler's Theorem, 

therefore 2 T^ = 2 j;=;?^j^+;>A+ ' ' - (2) 

dT dT 

Let us try to get -—^ and -^-^ indirectly. 

dq^ dp^ 

88 
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From Cl). !^^ = !^ + !i:^!2i + £S!2i+..., 

But from (2), 

2'^=p^'A+p,'A+..., (4) 

Subtracting (3) from (4), we get 

^' = -^'. (5) 

Again, we have from (1), 

"Pi hi ^Pi ^% "Pi 

From (2), 2^ = j,+^,|+^,|+ . . .. (7) 

Subtracting (6) from (7), we get 

dT 

^ = ?i- (8) 

The Lagrangian equation 

dt dqj^ dq^ 

dT 
becomes p^. -f — ^ = Q^.. (9) 

We have also 6^^ = ?^. (10) 

The equations of which (9) and (10) are the type are known 
as the Hamiltonian eqaations of motion. The so-called canonical 
form of the Hamiltonian equations is somewhat more compact 
and will be given later, in Chapter IV. 



^,^-^ = Q, 
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The 2n equations of which (9) and (10) are the type form 
a system ot 2n simultaneous differential equations of the first 
order, connecting the n coordinates j^, j'g' ' ' '» ?n' ^^^ ^^^ ^ com- 
ponent momenta p^, J^g' * * '' Pn^ ^^^ ^^^ time, and in order to 
solve for any one coordinate we must generally, as in the case 
of the Lagrangian equations, form and make use of the whole 
set of equations. 

In concrete problems there is usually no advantage in using 
the Hamiltonian forms, but in many theoretical investigations 
they are of importance. It may be noted that in the process of 
forming T^ from T^ , q^ is expressed in terms of the jt?'s and q% 
and thus equation (10) is anticipated. 

To familiarize the student with the actual working of the 
Hamiltonian forms, we shall apply them to a few problems 
which we have solved already by the Lagrange process. 

15. (a) The equations of motion in a plane in terms of 
polar coordinates (v. Art. 3, («)). 

Here ^^ = f [^ + ^0']- 

—i=:p=,mr, 
cr 



Whence r = £-% (1) 

m 

<^ = ^,' ' (2) 



and T = — \p* + ^\ 



dr mr^ 

dT 
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p^B<f> = 4>r&^. (4) 

Our Hamiltonian equations are 

A = r*. (8) 

If we eliminate p^ and j^^, we get 



?W 



[r - r20'] = R, 



our familiar equations. 

(J) Motion of a bead on a horizontal circular wire (v. Art. 5, (a)). 



Here T. = J «'^- 

— ^ = P# = matJ. 
cff 

mar 



T_ = 



^ P^' 



9 ' 2 
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Integrating, + -^=(7 = - 



P0 ma maV 

j__ 7» Vk 



ak Vkt -f wi 



(<?) The tractrix problem (v. Art. 6, (a)). 
Here i; = ^ [i* + a"^ - 2 a cos ^i^]. 

jt?^ = 771 [x — a cos 5(J], 
P0 = m \_a^d — a cos ^i]. 

Whence i = r-^- [jo^ + cos Op^'j^ (1) 

^ = — ^-^-f^ [« cos 6!p, + Pe] . (2) 

^P = H 2 • 2/1 [«X' 4-i?| -f 2 a cos ^i>^(,]. 

^ 2 wa* sm* ^ "- ^ ^ jtxt^j 

We get j9^ = .« sin ^, (3) 

Pe - ^^2 g^a [(«'i^x +i>.0 cos « + a (l-f cos^ e^^p^p,] = 0. (4) 

We have the condition 

x = nt (5) 

With (5), p„ = mln — acosdd^j 

» 
P0 = ma \a6 — n cos ^]. 

Substituting in (4), we get 

p^ — mna. sin Od = 0,, 

or mxi^O = 0. 
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Whence ^=C = -. 

a 

p^ = mn(l — cos tf), 
p^ = mn sin Od = sin tf = ^ sin ft 



a 



ii=^ 



a 
as in Art. 6. 

(rf) The problem of the two particles and the table with a 

hole in it (v. Art. 8, (<f)). 



m 



Here T.= ~[2i^ + (a-a:y^] 



We get 



and 



Whence 



if a; is small, as in Art. 8, (c?). 



Px- 


= 2 mi. 






Pi- 


-m(a — . 


xyd. 




• 


_ Px 
2m 






d: 


P» 
m(a — . 


-r 




T,- 


_ 1 \Pl 
2mV2 


' Pi 




\a-: 


zy 


« .1. 


Pi 


= mff, 




^' ' m(a-xy 






p» = 0. 






P» 


= C = ma ^ag. 




P.+ 


(a-xf 


g = mg. 






?r + ^-''^ 


x=0 






a 







\ 
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(«) The gyroscope (v. Art 10, (c)). 

^i = H'^(.^ + sin* Of*) + Cai cos <? + <^)']. 



dT. 






p^, = —^ = A sin» tf^ + C cos tf (^ cos ^ + <^), 
cf 

We get p^. = 0, (1) 

P* = 0, (2) 

P* ~ TTlTJa C(^* +-P*) cos - (1 + cos' 0) p^p^'] = mffa sin ^. (3) 
^ Sin u 

p^ = e'er, 

P^ ?_ [( cV + X") cos ^ - CLa (1 4- cos' ^)] = w^a sin tf, 

^v ex — Ca COS ^) ex COS ^ — Ca) . . ^ ,,. 

^^ = ^^^ , <^8/i ^ + ^^^ sin e. (4) 

^ sin* ^ "^ ^ ^ 

• • 

-^ cos + <j> = a, 

A sin' ^-^ + Ca cos ^ = X, 
as in Art. 10, (c?). 

16. The last two problems have a peculiarity that deserves 
closer examination. Let us consider Art. 15, (e). The kinetic 
energy in the Lagrangian form T^, and therefore* in the Ham- 
iltonian form T^, fails to contain the coordinates <f> and ^. 
Moreover, when either of these coordinates is varied, the 

* Since — ^ = (v. Art. 14), it follows that if a coordinate is missing in 

dQk dQk 

7^, it is missing also in Tp, 



or 
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impressed forces do no work. Hence two of our Hamiltonian 
equations assume the very simple forms 

which give immediately 

p^ = (7a, a constant, 

and p^ = L^ a constant. 

These enable us to eliminate p^ and p^ from a third Hamilto* 
nian equation (Art. 15, (^), (3)), which then contains only the 
third coordinate and its corresponding momentum pg. 

This same result might have been obtained just as well by 
replacing p^ and p^ in T^ by their constant values and then 
forming the Hamiltonian equations for in the regular way. 
So that if we are interested in only, and T^ has once been 
formed and simplified by the substitution of constants for p^ 
and p^y the coordinates <j> and yjr need play no further part in 
the solution. Should we care to get the values of these ignored 
coordinates, they can be found from the equations p^ = Car, 
p^ = Z, by the aid of the value of previously determined. 

In Art. 15, (cT), since p0 = O and pg = ma va^, the substitution 
of this value for pg in T^ enables us to solve the problem so 
far as a: is concerned without paying further attention to 0. 

To generalize, it is easily seen that if the Lagrangian form, 
and therefore the Hamiltonian form, of the kinetic energy fails 
to contain some of the coordinates of a moving system,* and 
if the impressed forces are such that when any one of these 
coordinates is varied no work is done, the momenta jt?^, p^^ * ' * 
corresponding to these coordinates are constant ; and that after 
substituting these constants for the momenta in question in the 
Hamiltonian form of the kinetic energy, the coordinates corre- 
sponding to them may be igno;:ed in forming and in solving 
the Hamiltonian equations for the remaining coordinates. 

* Coordinates that do not appear in the expression for the kinetic energy 
of a moving system are often called cyclic coordinates. 
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Unfortunately the ignored coordinates have to be used in 
forming T^, the Lagrangian form of the kinetic energy, and in 
deducing from it T^, the Hamiltonian form of the energy. 

Not infrequently this preliminary labor may be abridged 
considerably by using a modified form of the Lagrangian 
expression for the kinetic energy of the system, as we shall 
proceed to show. 

17. Routh's Modified Form of the Lagrangian Expression for 
the Kinetic Energy of a Moving System. In forming the 
Hamiltonian equations of motion (v. Art. 14) we first changed 
the form of T^ by replacing all the generalized velocities 
?i' ?2» • ' ' ^y their values in terms of the coordinates q^^ q^^ 

. . . and the generalized momenta p^^ p^, . . ., where p^ = — 2, 
i>, = T-~» etc. 

Let us now try the experiment of replacing in T^ one only 
of the velocities q^ by its value in terms of the corresponding 
momentum jp^, the coordinates q^j y^, • • ., and the remaining 
velocities q^<, q^^ • • •. 

Call T, thus changed in form, J^ . Of course 

T^ = T., and q,=^F(p^, q^, %^'-^ 9^ %^ ' ' 0- 

We have 'J^JJk + '^'A. 

dq^ dq^ dq^ 8q^ 

^ . dT. dT„ dq, a rm 

Transposmg, ^ = -^-.^^^^ = - [T,^_^^,J. 

Hence _ j^ = _[r,,-^jjj. 

is called the Lagrangian expression for the kinetic energy modified 
for the coch'dinate q^. 
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Our Lagrangian equation 

dt dq^ dq^ ^ 
becomes p^ — — ^ = Q^. (1 ) 

We have also o. = — — ^ • C2') 

It is noteworthy that (1) and (2) differ from the Hamil- 
tonian equations for q^ only in that the negative of the 
modified expression Mg^ appears in place of the Hamiltonian 
expression T^. 

Let us go on to the other coordinates. 

• 

^?, ^9. ^?i %, 
whence ^-!££i_« ?2x-l-rr -«al-^. 

89, 8q^ aj, aj/ 

The Lagrangian equation for q^ is therefore 

d dMg dMg 

and differs from the ordinary form of the Lagrangian equation 
only in that T^ is replaced by the modified expression M^^. 

In forming the modified expression it must be noted that q^ 
must be replaced by its value in terms of p^^ q^y ?8' * ' * ' ?i» 
9^, . . • , not only in T^ but in the term p^q^ as well. 

An advantage of the modified form is that when it has once 
been formed we can get by its aid Hamiltonian equations for 
one coordinate and Lagrangian equations for the others. 



48 MODIFIED LAGBAJNGIAN EXPRESSION [Art. 18 

The reasoning just given can be extended easily to the case 
where we wish Hamiltonian equations for more than one coordi- 
nate and Lagrangian equations for the rest. 

The results may be formulated as follows: Let Tp^,p^...,p^ be 
the form assumed by T^ when j^, ^j, • • • , 9^, are replaced by 
their values in terms of p^, P^y'-'y Pry ?r+i» ?r+2» • • •» in^ ?i» 
92» • • • » 9n' Then, if 

^\. 9,, •••.«.= ^Pv Pr '"^Pr- PAi - P2% Pr9ry 

we have equations of the type 

Pk r^^ — Vt» 

?* — '^ ' 

ifA<r + l; ^* 

if A > r. 

18. If we modify the Lagrangian expression for the kinetic 
energy for all the coordinates, 

^.,.",..= ^p -PA -M PrAn'y 

and we get Hamiltonian equations of the form 

.^^ ^^. (2) 

^Pk 

for all the coordinates, and as we have nowhere assumed in our 
reasoning that T^ is a homogeneous quadratic in the general- 
ized velocities, we can use these equations safely when the 
geometrical equations contain the time explicitly (v. Art. 7). 
If the geometrical equations do not involve the time, in which 
case T^ is a homogeneous quadratic in j^, Jj' * * * » 

2 ^i^Piii +M + • • • +Pn% 
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by Euler's Theorem ; and itf,,,-,,, = T^ - 2 T^ = - T^; and (1) 
and (2) assume the familiar forms 

A.+g=e. (3) 

dT 

It is important to note that the modified Lagrangian expression 
Mg^^ (?,, ■ •, g,. is not usually the kinetic energy of the system, although, 
as we shall see later, in some special problems it reduces to the 
kinetic energy. As we have just seen, when the time does 
not enter the geometrical equations, the completely modified 
Lagrangian expression (that is, the Lagrangian expression modi- 
fied for all the coordinates) is the negative of the energy. 

19. As an illustration of the employment of the Hamiltonian 
equations when the geometrical equations contain the time, let 
us take the tractrix problem of Art. 6, (a'). 

Here ^<? = ? \7^^ Jfo^^'-^ an cos dd^ 

and is not homogeneous in 6. 

P0 = —J = m \a^d — an cos ^], 

and & = -^ + -eoad. (1) 

am a 



^ 2 L a^mj 






-—2 = mn* sm cos ^ + - sm Op^ ; 
00 a 

n 

and we have p0 — mn^ sin cos sin 0p0 = ; (2) 

• a 

and (1) and (2) are our required Hamiltonian equations. Let 
us solve them. 
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From (1), p^^ma^d — mndoo^Oy 

whence p^ = ma^O + mna sin Od, 

Substituting in (2), 
ma^d+ mna sin 0d — mn^ sin cos tf — mna sin 5d + mn^ sin tf cos 0=Oy 
or d*=0, 

which agrees with the result of Art 6, (a'). 

EXAMPLES 

1. Work Art 6, (6'), by the Hamiltonian method. 

2. Work Exs. 1 and 2, Art. 6, by the Hamiltonian method. 

20. (a) As an example of the employment of the modified 
form, we shall take the tractrix problem of Art. 6 and modify 
for the coordinate x. 

We have (v. Art. 6) 

T^ = ^[d? + a^^^-2acos0x^']' 
p^ = m\x'-a cos 0&\ 

i = ^ + a cos 0d^ 
m 

2L TT? m ^'\ 

-77^ = 0, 

OX 

— ^ = ma' sin" 66 — a cos dp,, 

— -s = ma^ sin co» ^^ + a sin 6dp^. 
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We have for x the Hamiltonian equations 

p^ = Esin0, (1) 

i = ^ + acos0&, (2) 

m 

and for the Lagrangian equation 

ma' [sin^ 00 -h sin tf cos 0^- aco8 0p^ = O. (3) 

Of course (1), (2), and (3) must be solved as simultaneous 
equations, and we can simplify by the aid of the condition 

x=^nt, (4) 

Solving, we get ^ = 0, 

4 

^ = ^!L . (v. Art 6 and Art 15, (r)) 
a 

(6) As a second example we shall take the problem of the 
two particles and the table with a hole in it (v. Art. 8, (rf)) 
and modify for 0. 

We have T^ = 5[2i' + (a-a:)^]. 

p^ = —^ = mCa-' xYd, 

whence ^= ^'^^ ,, ■ (1) 

m(a — xy ^ ^ 



2L wi^Ca-a;)*] 



Our Hamiltonian equations are (1) and 

/, = 0. (2) 
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Our Lagrangian equation is 



m 



Ux-^ .^/' ^ = mg. (3) 



By (2), p^ = C=^maV^; 

whence (3) becomes 2 i H ^-—. = g^ r4) 

{a — xy 

as in Art. 8, (d). 

(c) As a third example we shall take the wedge and sphere 
problem of Ex. 3, Art 9, and modify for x. 

We have 

f^ = — -T = {M-\- m)x — my cos a. 

ex 

,^p^micosa^ 

M-^m ^ ^ 

^ jp^ - Tyiy cos^ g m a^ + ^ a 
"' 1{M^m) 2 a^ ^' 

^ p^-y^ycos^a m a^ + F ^ p^ -]- mpjf cos a 
2QM-\-m) 2 a^ ^ M-^m 

Our Hamiltonian equations are (1) and 

A = 0. (2) 

Our Lagrangian equation is 

[a^ -f ^ m^ cos^ al .. 7w»^ cos a . ,„^ 

w y ^-^ = mq sm a, (3) 

By (2), P.= C=Oi (4) 

, xoN i_ fa^-f-^^ m costal.. . .r^ 

whence (o) becomes ^ V = a' sm a, (5) 

as in Art. 9, Ex. 3. 
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(^d) As a fourth example we shall take the flexible-parallelo- 
gram problem (v. Art 9, (J)) and modify for <f>. 

We have ^^ = ? [2 J^<i>' -f «'^]. 

Our Hamiltonian equations are (1) and 

P* = 0. (2) 

Our Lagrangian equation is 

ma^0 = — mga sin 0^ (3) 

or ^"-f-^sin^ = 0. (4) 

a 

EXAMPLES 

1. Take the dumb-bell problem of Art. 8, (c), and modify for 0. 






mx = 0. 
my = 0. 
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2. Take the dumb-bell problem of Art 8, (c), and modify for 



^'=lk-^] 



^ ma^ff = 0. 
3. Take the gjrroscope problem, Art. 10, (c), and modify for 

^ -0 6=^^"^ C p^-^p^ COB 0)coB0 

^ _o, Jr=.P±:zP±^^. 

A0 — -^^^^^ — ^-^*^^^^ . .9^^ ^ = ^w^« sm ft 

^ sm^ 

21. We proceed to comment on jbhe problems of the pre- 
ceding section. 

(a) No one of them involves the tune explicitly in the 
geometrical equations, and therefore the kinetic energy T^ in 
all of them is a homogeneous quadratic in the generalized 
velocities. 

(6) The momenta, therefore, are homogeneous of the first 
degree in the velocities, and consequently the eliminated veloci- 
ties are homogeneous of the first degree in the corresponding 
momenta and the remaining velocities, and the energy ^^.j,,,... 
and the modified function J[f^^,^^,... are homogeneous quadratics 
in the introduced momenta and the velocities not eliminated. 
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(<?) In all the problems the coordinates for which we have 
modified the Lagrangian expression for the kinetic energy are 
cyclic. In all of them except the first no work is done when 
any one of the coordinates in question is varied. Consequently 
one of the Hamiltonian equations for that coordinate is of the 
form p^ = 0, and the momentum p^ = <?^ where ^^ is a constant. 
Therefore it is easy to express the energy ^p^, ,,.... and the 
modified expression Ji^^.,^,... in terms of the remaining coordi- 
nates, the corresponding velocities, and the constants c^^ <?,,•••, 
and when so expressed they are quadratics in the velocities but 
not necessarily homogeneous quadi*atics. When the modified 
function has been so expressed, it may be used in forming the 
Lagrangian equations for the remaining coordinates precisely 
as the Lagrangian expression for the kinetic energy is used, 
and the coordinates that have been eliminated may be ignored 
in the rest of the work of solving the problem unless we are 
interested in their values (v. Art. 16). 

(d) To generalize : If some of the coordinates of a moving sys- 
tem are cyclic^ and if the impressed forces are such that when any 
one of these coordinates is varied no work is done, the momenta 
corresponding to these coordinates are constant throughout the 
motion. The substitution of these constants for the momenta in 
the Lagrangian expression for the kinetic energy modified for the 
coordinates in question will reduce it to an explicit function of 
the remaining coordinates, the corresponding velocities, and the 
constants substituted, which will be a quadratic in the velocities 
but not necessarily a homogeneous quadratic. 

When the modified function has been so expressed, it may 
be used in forming the Lagrangian equations for the remaining 
coordinates precisely as the Lagrangian expression for the 
kinetic energy is used, and the coordinates that have been 
eliminated may be ignored in the rest of the work of solving 
the problem (v. Art. 16). 

In the important case where the system starts from rest, the 
constant momenta corresponding to the ignorable coordinates 
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being zero at the start are zero throughout the motion, and the 
modified expression is identical with the Lagrangian expression 
for the kinetic energy, which therefore is a function of the 
remaining coordinates and the corresponding velocities and is 
a homogeneous quadratic in the velocities (v. Art 24, (a)). 

(e) The fact that the Lagrangian expression for the kinetic 
energy modified for ignorable coordinates is expressible in terms 
of the remaining coordinates and the corresponding velocities 
and is a quadratic in terms of those velocities is often of great 
importance, as we shall see later. 

22. Let us take the gyroscope problem of Art. 10, ((?), and 
Art. 20, Ex. 8, and work it from the start, ignoring the cyclic 
coordinates ^ and '^. 

We have T.= i [^(^ + sin2^^')-h C(^costf •f<^y], and 
therefore <f> and -^ are cyclic coordinates. Moreover, no work 
is done when <f> is varied nor when -^ is varied, so that ^ and 
'^ are ignorable. 

b^ = 0, and p^ = 0, so that p^ = e^, and p^ = c^. 



We have p^ = —r=C(;fco80 + <^)=c^, \ 



p^ = —r = A sin^ 0yjr+C cos 0(;f cos 0-^^^=c^ 



u I ^1 C^2 — <^i COS 0) cos 
whence 6 = -i — -^-^ ^ . ^ / , 

^ C A 8111^0 

J : t?n — C. COS 

and ^ = ' .' 2 >. — 

A BUT 



T -^A^ I ^' I (^2-^1 COS g)n 

M -^A^^ ^^ (^2 - ^1 COS sy ] 
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Forming the Lagrangian equation for in the usual way, we 

^^""^ .s (^f + c|) cos e - c,c^ (1 + coa^ g) . . 

Ad - ^^^ ^-^ y ^ ^ = mija sni 6, 

A sm^ 

Ad . (^1 ~ ^2 cos ff) (e, — c, cos tf) . ^ 

or -4^ + -^^^^ , .^ ' — = mga sm ^, 

A sin*^ 

which is identical with (7), Art. 10, (c), 

EXAMPLE 

Work the problem of Art. 20, (J), ignoring the coordinate ft 

23. The problem of Art. 20, (rf), has an interesting peculiarity. 
Both coordinates are cyclic, and there is no term in the kinetic 
energy that is linear in <^, the velocity corresponding to the 
ignorable coordinate <^. Consequently the constant momentum 
j9^ is a constant multiple of <^, which is therefore itself a constant. 
This is true also of p^^ and of mF</)', the term in the energy 
which involves <^. T^^ and M^ must then differ from mo^^^ by con- 
stants, and as only the derivatives of M^ are used in forming the 
Lagrangian equation for ^, we have merely to disregard the terra 
mlc^^^ in the energy T^ and use what is left of T^ instead of M^. 

In cases like this we are able practically to ignore the con- 
tribution to the energy made by the ignored coordinate as well 
as the coordinate itseK, and of course we can conclude that the 
motion we may thus disregard has no effect on the motion we 
are studying, but that the two can go on together without 
interference. 

EXAMPLE 

Examine Exs. 1 and 2, Art. 20, from the point of view of 
the present article. 

24. (a) The wedge and sphere problem of Art. 20, (c), belongs 
to a very important class. Both x and y are cyclic coordinates, and 

X is ignorable. 

dT 

p^ = -~ = (Jtf 4- m^x — my cos a, p^ = 0. 

€/X 
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The momentum p^ is constant, and as it is initially zero 
(since the system starts from rest) it is zero throughout the 
motion, as is pjc. Consequently the kinetic energy T^ and the 
modified expression 31^ are identical, and as they are both 
homogeneous quadratics in y and p^ and do not contain y, they 
reduce to the form Ly^j where Z is a constant. 

Therefore our Lagrangian equation for y is Ly = mg sin a, 
and the sphere rolls down the wedge with constant acceleration. 
If we care only for the motion of the sphere on the wedge, we 
may then ignore x completely and yet know enough of the 
form of M^ to get valuable information as to the required 
^motion. Of course we know that the energy of the whole 
system can be expressed in terms of y, and if we are able by 
any means to so express it, we can solve completely for y with- 
out using the ignored coordinate x at any stage of the process. 

(6) As a striking example of ' this complete ignoration of 
coordinates, and of dealing with a moving system having an 
infinite number of degrees of freedom, let us take the motion 
of a homogeneous sphere under gravity in an infinite incom- 
pressible liquid, both sphere and liquid being initially at rest. 

From considerations of symmetry, the position of the sphere 
can be fixed by giving a single coordinate a;, the distance of 
the center of the sphere below a fixed level, and x is clearly a 
cyclic coordinate. 

The positions of the particles of the liquid can be given in 
terms of x and a sufficiently large number (practically infinite) 
of coordinates q^^ q^^ • • ., in a great variety of ways. Assume 
that a set has been chosen such that all the q^s are cyclic* 
Then, since gravity does no work unless the position of the 
sphere is varied, the ^''s are all ignorable. That is, for every one 
of them p^. = 0, and the momentum jt?^ = c?^, and since the system 
starts from rest the initial value of pj^ is zero, and therefore 
6'^ = 0. Jf,^,,^,..., the energy of the system modified for all the g^'s, 

* That this is possible will be shown later, in connection with the treatment 
of Impulsive Forces (v. Chap. Ill, Art. 36). 
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is then identical with the energy T^^.p,.... and must be expressi- 
ble in terms of the remaining coordinate x and the corresponding 
velocity x (v. Art. 21, (d)) ; and as a: is cyclic the energy of the 
system will then be of the form L^^ where Z is a constant. 
Forming the Lagrangian equation for a;, we have 

Lx = mg^ 

and we leam that the sphere will descend with constant 
acceleration. 

Of course this brief solution is incomplete, as it gives no infor- 
mation as to the motion of the particles of the liquid, and since 
we do not know the value of Z, we do not leam the magnitude 
of the acceleration. Still the solution is interesting and valuable. 

The energy of the moving liquid, calculated by the aid of 
hydromechanics, proves to be \rt\!^^ where rrJ is one half the 
mass of the Hquid displaced by the sphere (Lamb, Hydrome- 
chanics, Art. 91, (3)), and therefore the energy of the system 
is 1^ (m -j- m') i:*, and this agrees with our result 

25. Summary of Chapter II. The kinetic energy of a mov- 
ing system which has n degrees of freedom can be expressed 

in terms of the n coordinates q^^ ?2' ' * *' ?«' ^^^ ^'^ ^ general- 
ized momenta jt?^, jt?^, • • •, j9„, and when so expressed it is a 
homogeneous quadratic in the momenta if the geometrical 
equations do not involve the time (v. Art. 14), and is called 
the Hamiltonian expression for the kinetic energy. If :7^ is 
the Hamiltonian expression for the kinetic energy, and T^ the 
Lagrangian expression for the kinetic energy, 

-— -^ = — T-^ 5 and — ^ = Oj.. 

The work done by the effective forces in a hypothetical 
infinitesimal displacement of the system, due to an infinitesi- 

[dT'\ 

this be written equal to Qk^q^,^ the work done by the actual 
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forces in the displacement in question, there will result the 
differential Equation of the first order, 

The 2w equations of which this and 

given above, are the type are known as the llamiltonian equa- 
tions of motion for the system. 

If, in the expression Ti -p/i^ -p^q^ p^q^, q^, ^a^ • • > ?r» 

are replaced by their values in terms of p^, p^, • • ., jt?^, j^^^, 

•••» ?«' ?i» ?2' ' • '' ?»' *^® result Mg^ ^ is the Lagrangian 

expression for the kinetic energy modified for the coordinates 
?i' • • •» ?r' ^^^ ^^^® coordinates for which the expression has 
been modified, that is, for i<r-fl, we have Hamiltonian 

equations of the type 

. CM _ 

dq^ 

dM 
9k=-^ — 

For the remaining coordinates, that is, for k> r^ we have 
Lagrangian equations of the type 

d dM cM ^ , A ^ ^ 

rf"^^"a7 = ^*- (V. Art. 17) 

Whether we work from the Lagrangian expression T^ for the 
kinetic energy, or the Hamiltonian expression T^, or the modi- 
fied Lagrangian expression iM^^,..., we are in general led to a set 
of simultaneous differential equations whose number depends 
merely upon the number of -degrees of freedom in the moving 
system, and such that to solve one we must form and solve all. 

If, however, some of the coordinates are cyclic (v. Art. 16), 
and the impressed forces are such that when any one of them 
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is varied no work is done, the corresponding momenta are con- 
stant throughout the motion, and these coordinates are ignor- 
able in the sense that if constants are substituted for the 
corresponding momenta in Tp or -3lf,,...., the Hamiltonian equa- 
tions for the remaining coordinates in the former case (or 
the Lagrangian equations in the latter case) can be formed and, 
if capable of solution, can be solved without forming the equa- 
tions corresponding to the ignored coordinates. 

If the system starts from rest and there are ignorable coor- 
dinates, Jtf'g^,..., the Lagrangian expression modified for the 
ignorable coordinates, is identical with the kinetic energy of 
the system; and whether the system starts from rest or not, 
Mg^^... is a quadratic, but not necessarily a homogeneous quad- 
ratic, in the velocities corresponding to the coordinates which 
are not ignorable. 



CHAPTER III 

IMPULSIVE FORCES 

26. Virtual Moments. If a hypothetical infinitesimal displace- 
ment is given to a system, the product of any force by the 
distance its point of application is moved in the direction of 
the force is called the virtual moment of the force, and the sum 
of all the virtual moments is called the virtual moment of the 
set of forces. 

If the forces are finite forces, the virtual moment is the vir- 
tiuil work ; that is, the work which would be done by the forces 
in the assumed displacement. If the forces are impulsive forces, 
the virtual moment is not virtual work but has an interpreta- 
tion as virtual action, which we shall give later when we take 
up what is called the action of a moving system. 

27. For the motion of a particle under impulsive forces we 
have the familiar equations 

^ (^1 - ^0) = ^' 

are called the effective impulsive forces on the particle and are 
mechanically equivalent to the actual forces. 

If the point is given an infinitesimal displacement, 

is the virtual moment of the effective forces and of course is 
equal to the virtual moment of the actual forces. 

62 
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If the generalized coordinates of a moving system acted on 
by impulsive forces are y^, Jg' ' ' *' *^^ ^ displacement caused 
by varying q^ by hq^ is given to the system, 

where hqA represents the virtual moment of the effective forces. 

dx dx 



As in Art. 7, 



Hi Hx 



rrx. . ^x , dx m d ,^. 

1 hereiore mx -— = mx -7- = — — - ( ir ), 



and 8,^ = 



■-[CD.-©.: 



s?i- 



(where Pq^^q^ is the virtual moment of the impressed impulsive 
forces and Pq^ is called the component of impulse corresponding 
to q^ is our Lagrangian equation, and of course we have one 
such equation for every coordinate 5'^. 

Equation (1) can be written in the equivalent form 

28. Illustrative Examples, (a) A lamina of mass m rests on 
a smooth horizontal table and is acted on by an impulsive 
force of magnitude P in the plane of the lamina. Find the 
initial motion. 

Let (x^ y) be the center of gravity of the lamina, let 6 be the 
angle made with the axis of X by a perpendicular to the line 
of action of the force, and let a be the distance of the force 
from the center of gravity. 
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Then i'^ = I [i* + y« + i«^. (V. App. A, § 10) 

CPxX=<ipX=CPi)o = ^- 

If the axes are chosen so that x =0, y = 0, and d = 0, 
we have 

mx = 0, 



my-P, 




mk'd = aP. 




ties are 




m 


aP 



Hence the initial velocities are 

a; = 0, 

The velocities of a point on the axis of X at the distance b 
from the center of gravity are 

The point in question will have no initial velocity if 6 = . 

a 

It follows that the lamina begins to rotate about an instan- 
taneous center in the perpendicular from the center of gravity 

to the line of action of the force at a distance from 

a 

that line and situated on the same side of the line of the force 
as the center of gravity. This point is called the center of 
percussion. 
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(5) A wedge of angle a and mass M, smooth below and 
perfectly rough above, rests on a horizontal plane. 

A sphere of radius a and mass m is rotating with angular 
velocity ft about a horizontal axis parallel to the edge of the 
wedge and is placed gently on the wedge. Find the initial 
motion (v. Art. 9, Ex. 3). 

Take as coordinates Xj the distance of the edge of the wedge 
from a fixed axis parallel to it in the horizontal plane ; y, the 
distance of the point of contact of the sphere down the wedge ; 
and 0, the angle through which the sphere has rotated. 

We have 

both before and after the sphere is set down. After the sphere 

is set down, y — ad = 0. Before it is set down, x = y = and 

^ = ft. Since the sphere cannot slip, it exerts an impulsive 

force F mp the wedge, and an equal and opposite force F is 

exerted on it by the wedge at the instant the two bodies come 

in contact. 

dT- 
We have p^ = -^ = (^M+ m)x — my cos a:, 



ox 
dy 



Py = -^ = rn(^y-x cos a), 



Pg = — r = mf(^d. 
30 

Our equations are 

(^M-^- rn)x — my cos a = 0, (1) 

m(y — X cos a) = P, (2) 

7nA;'(^-ft) = -aP, ' (3) 

and we have also a6 = y. (4) 

From (2), (3), and (4), 

-—2— y - 2^ COS a = — - . (5) 

a a 
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From (1) and (5), 

QM-j-m) mco9 a \x = 11 cos a, 

and CM+ m) ^ m cos* «^ y = (-M'-f m) — ft. 

(tf) In Art. 8, (6), let the weight 4 m be jerked down with a 
velocity v. Find the initial motion. Take the coordinates x 
and y as in Art. 8, (6), and let P be the magnitude of the jerk. 





P. 


dx 


= m (8 a; - ^), 




Py 


dy 


= m (3 ^ — a;). 


Our equations 


are 


m(8 


x-y^=^P, 


and 




m(3 


y-x)=0. 


Whence 

• 






23 mi = 3 P, 
23 my - P. 


But 






3 P 

23 m 


and 






y = ^v. 



(d) Four equal rods freely jointed together in the form of a 
square are at rest on a horizontal table. A blow is struck at 
one corner in the direction of one of the sides. Compare the 
initial velocities of the middle points of the four rods. 

Let m be the mass and 2 a the length of a rod. Take a pair 
of rectangular axes in the table. Let (a:, «/) be the center of the 
figure at any time and 6 and <\> the angles made by two adjacent 
rods with the axis of X a:, «/, 0^ <^, are our generalized coordinates. 

The rectangular coordinates of the four middle points are 

obviously (x-aoo^<^,y-a sin <^), (1) 

(^x-\-a cos 6^ y -{-a sin ^), (2} 
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We have 



(x + a cos ^^ y + a sin <^), (3) 

(x — a cos 0, y — a sin 0). (4} 



Px = 4: mxy 

Let the values of x, y, ^, </>, be 0, 0, 0, — , before the blow is 
struck, and let P be the magnitude of the blow. 
Our equations are 4 wi = P, 

4 my = 0, 
2w(a'-f i^)^"=0, 
and 2m(a^ + A;^)<^ = aP. 

Whence y = 0, 

. P 



a; = 



im 

aP 



Let Vj, Vg, Vg, v^, be the required velocities of the four middle 
points. Then 

. . r 3a'-hA^ P 5P 

^1 ^^ 'P a^^yk^ 4m 8/w 

P 2P 

^2 = ^ = 1 = Q~' 

^ 4m 8 m 

J k^-a^ P IP 

P 2P 



V. = 2: = 



* 4m 8 m 



and Vj^ : ^2 : t^g : v^ = 5 : 2 : — 1 : 2. 
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29. General Theorems. Work done by an Impulse. If a particle 
(x^ y, z) initially at rest is displaced to the position (x -f Sa:, 
y '\-iy<, 2 -f S2), the displacement can be conceived of as brought 
about in the interval of time ht by imposing upon the particle 
a velocity whose components parallel to the axes are w^, v^, w^^ 
where hx = u^St^ hy = vJSt, and Sz = w^t 

If the particle is initially in motion with a velocity whose 
components are w, v, w^ the displacement in question could be 
brought about by imposing upon it an additional velocity whose 
components are obviously u^ — w, v^ — v, w^ — w. 

Let a moving system be acted on by a set of impulsive 
forces. Let m be the mass of any particle of the system ; 
P^, Pj,, P,, the components of the impulsive force acting on the 
particle ; w, v, w^ the components of the velocity of the particle 
before, and w^, v^, w^ after, the impulsive forces have acted. 

If any infinitesimal displacement is given to the system by 
which the coordinates of the particle are changed by Sj-, Sy, 
and S2, we have the virtual moment of the effective forces 
equal to the virtual moment of the actual forces (v. Art. 27); 
that is, 

2iw [(Wj — u)hx-\- (y^ — v)Sy'\- (^w^ — vf) &] 

= S [P> + P^hy + PM']' 

If the velocity that would have to be imposed upon the parti- 
cle w, were it at rest, to bring about its assumed displacement 
in the time ht has the components u^^ v^^ w^^ Sx = u^Sty Sy = v^Sty 
Sz = wJSt, and the equation above may be written 

2w [(Wj — w) Wg + (y^ — t^) Vj + (w^ — w") wj 

= 2[i*,P, + i;,P, + t^,P,]. (1) 

Interesting special cases of (1) are 

= I.[nP, + vP^ + wP,l (2) 
^^\u^P^^v^P^^w^p;\, (3) 
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the displacement used in (2) being what the system would 
have had in the time St had the initial motion continued, and 
that in (3) what it has hi the actual motion brought about by 
the impulsive forces. If we take half the sum of (2) and (3), 
we get 

or, a systerrCs gain in kinetic energy caused by the action of 
impulsive forces is the sum of the terms obtained by multiplying 
every force by half the sum of the initial and final velocities of 
its point of application^ both being resolved in the direction of 
the force. 

This sum is usually called the work done by the impulsive 
forces. 

30. Thomson's Theorem. If our system starts from rest, 
le = v = t^ = 0, and formulas (1) and (3), Art. 29, reduce re- 
spectively to 

Sw \u^u^ + v^v^ + w^w;\ = 2 [u^P^ 4- v^P^ -f w^P;\, (1) 

and 27W \ul + v^ -h w^-] = S \u^P^ + v^P^ + w^P;\. (2) 

But the first member of (1) is identically 

X I {["' + ''"- + <J + [«> + ^> + '"H 

and subtracting (2) from (1) we get 

X 1 {[< + *'' + <^ - [< + ^' .+ <] 

If ifcg, Vg, ^^2' ^^® ^'^^ components of velocity of the particle m in 
any conceivable motion of the system which could give the 
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points of application of the impulsive forces the same velocities 
that they have in the actual motion, then the second member of 
(3) is zero and we get Thomson^s Theorem: 

If a %y%tem at rest is set in motion by impulsive forces^ its kinetic 
energy is less than in any other motion where the velocities of the 
points of application of the forces in question are the same as in the 
actual motion^ by an amount equal to the energy the system would 
have in the motion, which, compounded with the actual motion, would 
produce the hypothetical motion. 

31. Bertrand's Theorem. If Q^, Q^, and Q^ are the com- 
ponents of the impulsive force which would have to act on the 
particle 7n of the system considered in Art. 29 to change its 
component velocities from w, v, w^ to w^, v^, w^^ formula (3), 
Art. 29, gives us 

Stw [(Wjj - w) Wj -f (Vjj - v) Vj + 0^2 - w') w^ 

^-^Iu^Q^ + v^Q^-^w^Q:]. (1) 

Subtracting (1) from (1) in Art 29, we get 

2m l(u^u^ + v^v^ + w^w;) - (ill + «^2 + «^2)] 

= SK(P, - gj + v^(I\ - Q,) + w^iP^ - Q,)]. (2) 

The first member of (2) is (v. Art. 30) identically 

X ? {[< + ^' + '"?] - ["' + < + <'^ 

If the second member of (2) is zero, as will be the case if the 
Q-forces differ from the P-forces only by the impulsive actions 
and reactions due to the introduction of additional constraints 
which have no virtual moment in the hypothetical motion into 
the original system, we have Bertrand's Theorem: 

If a system in motion is acted on by impulsive forces^ the kinetic 
energy of the subsequent motion is greater than it would be if the 
system were subjected to any additional constraints and acted on 
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by the same impuUive forces^ by an amount equal to the eneryy it 
would have in the motion whichy compounded tvith the first motion^ 
would give the second.* 

32. By the aid of Thomson's Theorem many problems in- 
volving impulsive forces can be treated as simple questions in 
maxima and minima. 

(a) If, for example, in Art. 28, (a), instead of giving the 
force P we give the velocity v of the foot of the perpendicular 
from the center of gravity upon the line of the force, so that 
y -\-ad = Vy then to find the motion we have only to make the 

energy T^ = — [a? -\-y^ + 1^6^'\ a minimum. 

a 

dT- 

ex 

dT' }? 



dy a 

Whence i; = 0, 



W 



y^^f^'' 



^= ^ 



a^-vie 



V, 



and these results agree entirely with the results obtained in 
Art. 28, (a). 

* Gauss s Principle of Least Constraint : If a constrained system is acted on 
by impulsive forces, Gauss takes as the measure of the "constraint" what is 
practically the kinetic energy of the motion which, combined with the motion 
that the system would take if all the constraints were removed, would give the 
actual motion. 

It follows easily from Bertrand's Theorem that this "constraint" is less 
than in any hypothetical motion brought about by introducing additional con- 
straining forces (v. Routh, Elementary Rigid Dynamics, §§ 391-893). 
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(6) In Art 28, ((?), since x = r, the energy 

To make Tk a minimum, we have 



dTi, 


w[3y- 


t.]= 


0, 


y^ 


V 

V 







and the problem is solved. 

{<•) In Art 28, (rf}, let t)j = i + a^ be given. 

Tlieu ^i = f [-1 (", - ««^)* + 4 y* + 2 (a» + F)(<^ + <^')]. 

^ = 4 wy = 0, 

ad 

^' = »i[2(a» + F)<^ - 4 a(«)j - a<^)] = 0. 
dd) 

y = o, 
6 = 0, 

I 2av^ 

9 = 



x = 



a' + A" 2 

-—-,^ = ,^ = ,^ = -,^, 

. , -a' + F 1 

and ,, = a:-«^ = ^^j— ^t.,=--t.. 

33. In using Thomson's Theorem we may employ any vahd 
form in the expression for the energy communicated by the 
impulsive forces. For instance, in the case of any rigid body, 

is permissible and is much simpler than the corresponding form 
in terms of Euler's coordinates. 
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Take, for example, the following problem: An elliptic disk 
is at rest. Suddenly one extremity of the major axis and one 
extremity of the minor axis are made to move with velocities 
U and V perpendicular to the plane of the disk. Find the 
motion of the disk. 

Let us take the major axis as the axis of X and the minor 
axis as the axis of Y. 

We have, then, T=-^[f + f + e + Aa>l-\- B<ol+ Co)']. 

By the conditions of the problem, since the components of the 
velocity of the point (a, 0, 0) are 0, 0, Z7, and those of the 
point (0, by 0) are 0, 0, F, we have 

i:=0, 

^ + awg = 0, 

i — a(D^ = U\ 

and X — h(o^ = 0, 

y = 0, 

Hence r = |p + ^(F-i/ + |(Cr_iyj, 



or 



T ='^\^ + \iV - zf + \iu - zy'^. 



mh 1 ^, ma 

smce A = — — and B = — -- • 

4 4 

z = \(u+vy 

We have also ®i = ^^ (5 ^— ^)» 

bo 
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EXAMPLE 

One extremity of a side of a square lamina is suddenly made 
to move peq)endieular to the plane of the lamina with velocity 
U^ while the other extremity is made to move in the plane of 
the lamina and perpendicular to the side with velocity F. 

Show that the center will move with velocity — perpendicular 

to the plane, and with velocity — v 2 in the plane, toward the 
corner on which the velocity V was impressed. 

34. If our system does not start from rest, it is often easy 
to frame and to solve a problem in which the system is mitially 
at rest and is acted on by the same impulsive forces as in the 
actual problem, and where consequently the resulting motion, 
compounded with the actual initial motion, will give the actual 
final motion. 

For example, consider the following problem: A sphere 
rotating about any axis is gently placed on a perfectly rough 
horizontal plane. Find the initial motion. Here, in the actual 
case, the lowest point of the sphere is immediately reduced 
to rest. 

Take rectangular axes of X and F, parallel to the plane and 
through the center of the sphere. Let Xl^, Xl^, ft^, be the com- 
ponent angular velocities before, and o)^, cd^^, cd^ after, the sphere 
is placed on the plane. Let i, y, be the velocities of the center 
of the sphere. 

Then, in the actual case, x — aWy = and y + cuo^ — are our 
given conditions. Initially the velocities of the lowest point of 
the sphere are — aft^ and aO^. If the sphere were at rest, the 
impulsive force which in the actual case destroys these velocities 
would give to the lowest point the negatives of these velocities ; 
that is, ally and — aO^. Let us then solve the following 
auxiliary problem: A sphere is at rest. Suddenly the lowest 
point is made to move with velocities aO,^ and — all^, parallel to 
a pair of horizontal axes. Find the initial motion of the sphere. 
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If u and V are, respectively, the x component and the y 
component of the velocity of the center, and g)^, g)^, g)^, are the 
angular velocities. 

The velocities of the lowest point are w — aw^, v-f aa)^, but 
they were given as aft^ and — aO,^. 

Therefore u — aco^ = aH^, 

r = I K [K + ^.y + ('"i + "x)' + *" K + «l + «.*)]}• 

To make this a minimum, we have 

dT 



da) 

1 



= m [a\co^ + ft,) + A^G)J = 0, 






dT 

- — = TWJfe^ffl = 0. 

3(0 « 

Hence « = - ,^^ ft,, 

^ rt -I- At 

<»- = ^ IT, 11. 



2 a^ + F ^' 



0) =0. 

3 



Compounding these with the initial angular velocities in the 
actual problem, ft^, ft^, ft^, we get 



"^x- 


a^ + e "'' 


Q) — 

y 


a' + ie "' 



ft). 



= ft,. 



These equations, together with x — atOy = and y + og)^ = 0, 
completely solve the original problem. 
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35. A Problem in Fluid Motion. Let us now consider an inter- 
esting application of tlie principles of this chapter which was 
made by Lord Kelvin to a problem in fluid motion. 

It is shown m treatises on hydromechanics that if an, incom- 
pressible, frictionless, homogeneous liquid, either infinite in 
extent or bounded by any finite closed surfaces fixed or moving, 
and with any rigid or flexible bodies immersed in it, is moving 
imder the action of conservative forces (v. Chap. IV) and has 
ever been at rest, the motion will be what is called irrotational. 
That is, if a-, y, 2, are the rectangular coordinates of any fixed 
point in the space occupied by the liquid, there will be a func- 
tion <^(a;, y, z) such that if i^ t;, t^, are the components of the 
velocity of the liquid at the point 2:, y, 2, 

dS d(b cd) 

OX cy cz 

The function <^ is called the velocity-potential function. 

Since throughout the motion the liquid is always supposed 
to be an incompressible continuum^ v^ v, w^ must satisfy the equa- 
tion of continuity for an incompressible liquid, 

dx dy dz 
and therefore <^ satisfies Laplace's equation, 



^.?!^ 4.^ = 

da^ "^ dy^ ^ dz' ' 



and will be uniquely determined except for an arbitrary constant 

term if the value of --^, the velocity normal to the surface, is 

dn 

given at every point of the boundary of the liquid, however 
irregular that boundary. Therefore the actual motion at every 
point of the liquid at any instant is uniquely determined, if the 
motion is irrotational, when the normal velocities at all points 
of the boundary are given. 

We wish now to prove that the kinetic energy of the actual 
motion is less than that of any other motion, not necessarily 
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irrotational, consistent with the equation of continuity and with 
the actual normal velocities at the boundary. 

If tiy V, w, are the velocities at (xy y, z) in the actual motion, 
let u-\-aj v-f /8, t^-f 7, be the velocities in the hypothetical 
motion, and let v^ be the actual normal velocity at any point 

of the boundary. Then we have ■^ — lu-j-mv-\-nw = v,, where 

on 

Z, w, w, are the direction cosines of the normal, and w, v, Wj are 
the components of the velocity at the point in question. 
In the hypothetical motion, 

?(w + a) + 7W (t; + /8) -f- n (w -f 7) = v^ 

at the same point; therefore 

la'\-m^-hny = (1) 

at every point of the boundary. 

As the hypothetical velocities as well as the actual velocities 
must obey the law of continuity, 

dCu + a) ^ d(y-h^} ^ ^(^^ + 7) ^0 
dx dy dz 



^ 1 f ca dS dy ^ 

and therefore -- -f -^ + -i = 

ox oy 02 

at every point in the bounded space. 
If T is the energy of the actual motion, 



(2) 



-i///t 



w^ -h v^ + w^^ dxdydzy 

where p is the density of the liquid, and where the volume 
integral is taken throughout the' space filled by the liquid. 
If jT' is the energy of the hypothetical motion, 

T^ = ^ /YTEC^ + «)' + (^ + ^y + (^ + 7)'] dxdy dz 
= T ^^ CCCla' + ^ + r/]dxdydz 

+ P I j j [«^^ + /3v + ywl dxdydz. 
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By the aid of Green's Theorem we can prove that 

[aw H- /8v -f yw'] dx dy dz = 0. 



Ml 

We have III -r-dxdi/dz= I UcosadS 



(v. Pekee, Newtonian Potential Function, p. 92 (143)), where 
the volume integral is taken throughout any bounded space 
and the surface integral over the boundary of the space, cos a 
being the x direction cosine of the normal to the boundary. 

Now au = a-;^ = — ('a<p)— 07-, 

ex ex ex 

and lit audxdydz = j la^dS — III <f> — dxdydz. 

In like manner, 

III ^vdxdydz= I mfi^dS — III <^ — dxdydz^ 

and III 'iw dx dydz= I ny^dS — I 1 1 <f> -^ dx dy dz. 

Hence III [au -f ^y -|- ytv^ dx dy ds 

But the surface integral vanishes by (1), and the volume 
integral vanishes by (2). Therefore the energy T' is greater 
than the actual energy T. 

It follows that the irrotational motion of any frictionless 
incompressible homogeneous liquid under the action of con- 
servative forces is at every instant identical with the motion 
which would have been suddenly generated from rest by a set of 
impulsive forces applied at points in the boundary of the liquid 
and such that they would suddenly give all the points of the 
boundary the normal velocities that these points actually have 
at the instant in question (v. Thomson's Theorem in Art. 30). 



^ + ^ + ^Ur%&. 
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36. If now we have a liquid contained in a material vessel 
and containing immersed bodies, a set of generalized coordinates 
?i' ?2' ' • *' ?«' ^^^ ^ chosen, equal to the number of degrees of 
freedom of the material system formed by the vessel and the ' 
immersed bodies, and the normal velocity of every point of the 
surface of the vessel and of the surfaces of the immersed bodies 
can be expressed in terms of the coordinates q^^ q,^^ ••'»?»» ^^^ 
the corresponding generalized velocities j^, q^^ • • •, q„. 

We can now choose other ihdependent coordinates q{^ Ja, • • •, 
practically infinite in number, which, together with our coordi- 
nates j'j, ?2, • • •, g'„, will give the positions of all the particles 
of the liquid. 

Suppose the system (vessel, immersed bodies, and liquid) at 
rest. Apply any set of impulsive forces, not greater in num- 
ber than 71, at points in the surface of vessel and of immersed 
bodies and consider the equations of motion. For any of our 
coordinates qj^ we have the equation 

where p'^. is the generalized momentum corresponding to 9^., 
since (as in varying q'^ no one of the coordinates 9^, Jg' * * *' 5'n' 
is changed, and therefore no one of the impulsive forces has its 
point of application moved) the virtual moment of the com- 
ponent impulse corresponding to q^ is zero. 

As the actual motion at every instant could have been gen- 
erated suddenly from rest by such impulsive forces as we have 
just considered, the momentum pl^ is zero throughout the actual 
motion ; and the impressed forces being by hypothesis conserv- 
ative, and the liquid always forming a continuum, no work is 
done when 9^ is varied. Consequently, in the Hamiltonian 

dT dT 

equation »^ -f- —-7 = ^^, p'i, = Q and ^^ = 0, and therefore -—7 = 0. 

^k Hk 

Hence every coordinate q^^ is cyclic, and it is also completely 

ignorable. The energy modified for the coordinates q'j^ is identi- 
cal with the energy, which, being free from the coordinates q'j^ 
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and the moment^/?^, is expressible in terms of the n coordinates 

'/i» %"* ' • *' ^A» ^^^ ^^^® corresponding velocities q^^ q^j • • •, j„, 
and is a homogeneous quadratic in terms of these velocities 
(v. Art. 24, («)). 

37. Summary of Chapter III. In dealing with problems in 
which a moving system is supposed to be acted on by impul- 
sive forces, we care only for the state of motion brought about 
by the forces in question, since on the usual assumption that 
there is no change in configuration durmg the action of the 
impulsive forces, we are not concerned with the values of the 
coordinates but merely with the values of their time derivatives. 

The virtual moment (v. Art. 26) of the effective impulsive 
forces in a hypothetical infinitesimal displacement of the system, 
due to an infinitesimal change Sq^^ in a single coordinate q^j is 

If either of these is written equal to -Pj.S^'^, tbe virtual moment 
of the actual impulsive forces in the displacement in question, 
we have one of the equivalent equations 






The n equations of which either of these is the type are n 
simultaneous linear equations in the n final velocities (5'j)i, 

(5^2) 1' ' * *' ^^^ ^ ^^® configuration and the initial state of 
motion are supposed to be given, no integration is required, 
and the problem becomes one in elementary algebra. 

A skillful use of Thomson's or of Bertrand's Theorem re- 
duces many problems in motion under impulsive forces to 
simple problems in maxima and mimima. 



CHAPTER IV 

CONSERVATIVE FORCES 

38. If X, F, Z, are the components of the forces acting 
on a moving particle (coordinates x^ y, 2), the work, TT, done 
by the forces while the particle moves from a given position 
^0' (^0' ^0' ^0)' ^^ ^ second position P^, (a;^, y^, 2^), is equal to 






p. 
{Xdx + rrf^ + Zdz\ ; 



and since every one of the quantities X, F, and Z is in the 
general case a function of the three variables x^ y, 2, we need 
to know the path followed by the moving particle, in order to 
find W, Let f{x^ y, 2) = 0, <^ (a;, y, 2) = 0, be the equations of 
the path. We can eliminate z between these equations and then 
express y exphcitly in terms of a;; we can then eliminate y 
between the same two equations and express z in terms of x\ 
and we can substitute these values for y and 2 in X, which will 

I Xdx can be 

found by a simple quadrature. By proceeding in the same way 

with Y and Z we can find I Ydy and I Zdz^ and the sum 

of these three integrals will be the work required. 

It may happen, however, that Xdx -\- Ydy + Zdz is what is 
called an exact differential^ that is, that there is a function 
U=(t>(Xy y, 2) such that 

I^.X, ^^=Y, ^^=z. 

OX cy oz 

81 
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Since the complete differential of this function is 

ex cy cz 

or Xdx -h Ydy + Zdz^ we have 

[Xdx + Ydy + Zdz'] = <^ (a:, y, 2) = ?7, 



/I 



and r \xdx+ Ydy -i- Zdz'] 

and in obtaining this result we have made no use of the path 
followed by the moving particle. 

When the forces are such that the function U=<\>(x^ y^ z^ 
exists, they are said to be conservative^ and U is called the 
force function. 

We can infer, then, that the work done by conservative 
forces on a particle moving by any path from a given initial 
position to a given final position is independent of the path and 
is equal to the value of the force function in the final position 
minus its value in the initial position. 

If instead of a moving particle we have a system of particles, 
the reasoning given above applies. 

Let (r^, y^., Zf!) be any particle of the system, and JE^., Yj^^ Z^^ 
be the forces applied at the particle. Then the whole work, TF, 
done on the system as it moves from one configuration to 
another is equal to 

If there is a function TJ ^ ^{p^^i ^^g, • • •, a;;,.' ' * '» V-^ ^2' * * *' 
y^t? • • •<) ^1? 25^, • • ., 2^, • • •) such that 

^-]r ^_T- ^-y 

then 2 \_Xj^ dxj^ + Yj^ dyj^ + Zj^ dzj^ is an exact differential and 

U= 4>(x^i ^Jg, • • •, y^, ^2, • • •, ^j, ^2' • * •) ^'^ ^^® indefinite integral 
and is the force function ; the forces are a conservative set ; and 
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the work done by the forces as the system moves from a given 
initial to a given final configuration is equal to the value of 
the force function in the final configuration minus its value in 
the initial configuration, no matter by what paths the particles 
may have moved from their initial to their final positions. 

It is well known and can be shown without difficulty that 
such forces as gravity, the attraction of gravitation, any mutual 
attraction or repulsion between particles of a system which for 
every pair of particles acts in the line joining the particles and 
is a function of their distance apart, are conservative ; while such 
forces as friction, or the resistance of the air or of a liquid to 
the motion of a set of particles, are not conservative. 

The negative of the force function of a system moving under 
conservative forces is called the potential energy of the system, 
and we shall represent it by V, 

If we are dealing with motion under conservative forces and 
are using generalized coordinates, and the geometrical equations 
do not contain the time, we can replace the rectangular coordi- 
nates of the separate particles of the system in the force function 
or in the potential energy by their values in terms of the gener- 
alized coordinates g'j, 5^3' ' ' *' ?»' ^^^ ^^ ^^^ \h\i^ get U^ and con- 
sequently F, expressed in terms of the generalized coordinates. 

If U is thus expressed, Q^^qj, (the work done by the impressed 
forces when the -system is displaced by changing qj^ by hq^^ is h^JJ 

and therefore is approximately - — hq^^ or — •;;— hq^^ and hence 

dq^ cqj^ 



«.= 



dq^ "" dqj^ 



39. The Lagrangian and the Hamiltonian Functions. If the 
forces are conservative, our Lagrangian equation 

dt dq, dq, ^' ^ > 

may be written — -^ = -_i - ^ , (2) 

dt dq^ dqj^ dq^^ 
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where V is the potential energy expressed in terms of the coordi- 
nates j'j, ^ji • • '1 and not containing the velocities j^, jg, • • •• 

If L=T^- r, (3) 

£ is an explicit function of the coordinates and the velocities 
and is called the Livjranyian function. 

Obviously ^- = ^, and ;;— = —?-- — 

cq^ cq^ cq^ dqj^ dq^ 

Hence our Lagrangian equation (1) can be written very neatly 

dt\cqj dq. 

If our forces are conservative, and we are using the Hamil- 
tonian equations and express the kinetic energy T^ in terms of 
the coordinates and the corresponding momenta, and if we let 

If=T^-hV. (5) 

I£ is called the Hamiltonian function and is a function of q^^ q^^ 

..., q^, and ;?^, ;>2' • • •' A- 
Our Hamiltonian equations 



Pk^ 



p 



can now be written 



dq^ 


- Vt. 


9k 




^Vi. 


dH 


dt 


hk 


dt 


dH 



(6) 



(7) 

and these are known as the Hamiltonian canonical equations. 

If our forces are conservative, and we are using instead of 
the kinetic energy the Lagrangian expression for the energy 
modified for some of the coordinates q^^ q^^ • • •» Jr (y* ^^' 17)» 
and if we let 4> = .!/„.... .,,- F, (8) 
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we have for any coordinate q^^j^ the Lagrangian equation 



-- - 



and for any coordinate q^^^. the pair of Hamiltonian equations 

40. The Lagrangian function L is the difference between the 
kinetic energy T^ (expressed in terms of the coordinates j^, q^^ 
• . •, 5'^, and the velocities q^, q^^ • • •, q^, and homogeneous of the 
second degree in terms of the velocities) and the potential 
energy V (expressed in terms of the coordinates j^, q^, • . ., q^^. 

The Hamiltonian function If is the sum of the kinetic 
energy T^ (expressed in terms of the coordinates j^, q^y • . ., q^, 
and the corresponding momenta jt?^, jo^, • . ., jo„, and homogen- 
eous of the second degree in terms of the momenta) and 
the potential energy V (expressed in terms of the coordinates 

The sum of the kinetic energy and the potential energy, 
however expressed, is sometimes called the total energy of the 
system, and we shall represent it by E^ so that 

^=T+F. (1) 

The function <I> of the preceding section is the difference 
between the kinetic energy (expressed in terms of the coordi- 
nates q^^ jg, • • • , q^^ the momenta jt?^, p^, • • • , jt?^, and the 

velocities J^+i? ?r+2» • • '» ?«)? niinus the terms 

(similarly expressed), and the potential energy (expressed in 
terms of the coordinates q^^ q^^ • • • , q^. We shall call it the 
modified Lagrangian function. 
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It is to be observed that all the terms of <I> except those 
contributed by the potential energy are homogeneous of the 
second degree in the momenta introduced and the velocities not 
eliminated by the modification. 

41. In dealing with the motion of a system under conserva- 
tive forces, we may form the differential equations of motion 
in any one of three ways, and the equations in question are 
practically given by giving a single function — X, the Lagrangian 
function, or H^ the Hamiltonian function, or <I>, the modified 
Lagrangian function.* 

Every one of these functions consists of two very different 
parts : one, the potential energy Fi which depends merely on the 
forces, which in turn depend solely on the configuration of the 
system; the other, the kinetic energy T or the modified 
Lagrangian expression M^ ..., either of which involves the veloci- 
ties or the momenta of the system as well as its configuration. 

If we are using as many independent coordinates as there are 
degi'ees of freedom, a mere inspection of the given function, 
will enable us to distinguish between the two functions of 
which it is formed, the potential energy or its negative being 
composed of all the terms not involving velocities or momenta. 

If, however, we are ignoring some of the coordinates (v. Arts. 
16, 21, and 24) and are using /T (the Hamiltonian function) or ^ 
(the Lagrangian function modified for the ignored coordinates), 
tlie portion contributed to H by T^ or to <I> by J!f (the modified 
expression for the kinetic energy) is no longer necessarily a 
homogeneous quadratic in the velocities and momenta (v. 
Art. 21) and may contain temis involving merely the coordi- 
nates and therefore indistinguishable from terms belonging to 
the potential energy ; and consequently the part of the motion 
not ignored would be identical with that which would be pro- 
duced by a set of forces quite different fronj the actual forces. 

* Indeed, for equations of the Lagrangian type, any constant multiple of L 
or 4» will serve as well as L or 4>. 
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We may note that the last paragraph does not apply if the 
system starts from rest, so that the ignored momenta are zero 
throughout the motion (v. Art. 24, (a)). 

Let us consider the problem of Art. 8, (c?), where the potential 
energy V is easily seen to be — mgx. 

If we are using the Lagrangian method, as in Art. 8,<rf), 

we have _ 

^ = f [2^ + («-^)^ + 2^r]. (1) 

If we use the Hamiltonian method, as in Art. 15, (cf), we 
have 1 [v^ t?2 -I 

If we use the Lagrangian method modified for ^, as in 
Art. 20, (6), we have 

^^^hd^ P^' ,, + 2. J. (3) 

2L nr(a^xy J 

A mere inspection of any one of these three functions enables 
us to pick out the potential energy as F = — niffx. 

If, however, we are ignoring ^, as in ^Vrt. 22, Example, we 

z\_ m\a — xy J 

and here, so far as the function <I> shows, the potential energy 



e" 



may be — mgx^ as, in fact, it is, or it may be —^ — mgx. 

A Vfi \0l/ — Xj 

Indeed, the hanging particle moves precisely as if its mass were 
2 m and it were acted on by ^ force having the force function 

^= - — -2 + mgxi 

I m {a — xy 

that is, a force vertically downward equal to mg ^ ^_- 

m^a — xy 

If, as in many important problems (v. Chap. V), we are unable 

to discern and measure the impressed forces directly and are 

attempting to deduce them from observations on the behavior of 



(j^ 
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a complicated system, wliich for aught we know may contain 
undetected moving masses, the fact that we cannot discriminate 
with certainty between the terms contributed to the modified 
Lagrangian function by the kinetic energy of the system and 
those contributed by the potential energy may lead to entirely 
different and equally plausible explanations of the observed 
phenomena (v. Art. 51). 

42. Conservation of Energy. If we are dealing with a system 
moving under conservative forces, the coordinates q^^ q^, • • •, q^^ 
are functions of t, the time, as are also the velocities q^^ ^a' * ' *' ?»• 
Therefore Fi the potential energy, and T, the kinetic energy, are 
functions of the time, as is their sum, the Hamiltonian function H. 

Let us find -— • 

dt 

As ^depends explicitly on the coordinates q^, ?2' ' * "' ?«' ^^^ 
the momenta p^, J^j' * * *' i^«' 

dt ^[Sq^ dt dp,, dt] ^ ^ 

But by our Hamiltonian equations (v. Art. 39, (6) and (7)), 

^ = _?^, and %=.^, 
dt dqj^ dt dp^ 

and (1) becomes 

dt ^ [dq, dp, dp, dq,\ • ^ ^ 

Therefore T-\-V = ff==h, (3) 

where A is a constant. 

Hence in any system moving under conservative forces the 
sum of the kinetic energy and the potential energy is constant 
during the motion. 

This is called the Principle of the Conservation of Energy. 

Since by (3) any loss in potential energy during the motion 
is just balanced by an increase in the kinetic energy, and the 
loss in potential energy is equal to the work done by the actual 
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forces during the motion, our principle is a narrower statement 
of the familiar principle : If a system is moving under any forces^ 
c<mservative or not, the gain in kinetic energy is always equal to 
the work done by the actual forces. 

43. Hamilton's Principle. Let a system move under conserv- 
ative forces from its configuration at the time t^ to its con- 
figuration at the time t^. We have 

d dL dL .^ 

dt C(i^ cq^ 

where X, the Lagrangian function, is equal to T—V. 

Suppose that the system had been made to move from the 
first to the second configuration so that the particles traced 
slightly different paths with slightly different velocities, but so 
that at any time t every coordinate qj^ differed from its value 
in the actual motion by an infinitesimal amount, and so that 
every velocity q^^ differed from its value in the actual motion 
by an infinitesimal amount, or (using the notation of the 
calculus of variations)* so that hqj^ and hq^, were infinitesimal ; 
and suppose it had reached the second configuration at the 
time t^. Then, if at the time t the difference between the value 
of X in the hypothetical motion and its value in the actual 
motion is SX, ^„ _ ^ 

Now, at the time f, , 






dL^. dL d 



dt\cq, ^7 dtdq, ^* 



* For a brief introduction to the calculus of variations, see Appendix B. 
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Therefore, by (2), «Z = | J) [?^ «?*]' 

and fydt = Bj^'Ldt=[x'4M': 

Since the terminal configurations are the same in the actual 
motion and in the hjrpothetical motion and the time of transit 
is the samcy Sq^ = when t = t^ and when t = t^, and 

S f'Ldt=S f\T-V^dt = 0. (3) 

But (3) is the necessary condition that L should be either a 
maximum or a minimum and is sometimes stated as follows : 
When a system is moving under conservative forces^ the time inte- 
gral of the difference between the kinetic energy and the potential 
energy of the system is " stationary " in the actual motion. This 
is known as Hamilton's principle,* 

44. The Principle of Least Action. If the limitation in the pre- 
ceding section that '* in the actual and the hypothetical motions 
the time of transit from the first to the second configuration is 

* Hamilton's principle plays so important a part in mechanics and physics 
that it seems worth while to obtain a formula for it which is not restricted to 
conservative systems. 

We shall use rectangular coordinates, and we shall make the hypothesis as to 
the actual motion and the hypothetical motion which has been employed above. 

For every particle of the system we have the familiar equations 

mi = X, my = F, TTw; = Z. 
Since r = ^ I (*^ + 1^2 + i^), 

dt at at J 

But ?7ix — 5x = m — (xto) — mxbx = — (mxbx) — Xbx, 
dt Ot at 

Hence 5T + 2 [Xdx + Ydy + Z8z] = - 2m [x8x + ^5y + iSz], 
and r '{5T + 2 [XSx + Y^ + ZSz^jdt = 2 [mxdx + i^V + iSz]'.^ = 0. 

If the system is conservative, 2[X5x + Ydy + Zdz] = 5Z7 = — 5T, and we 
get the formula (3) in the text. 
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the same " be removed and the variations be taken not at the 
same time but at arbitrarily corresponding times, t will no longer 
be the independent variable but will be regarded as depending 
upon some independent parameter r. We now have (v. Art 43) 

dL ^. dL d ^ . dL d ^. 
and ^^ = |S[^^*]-|^'-2;[&g] 



smce 



dL dT^ 



Hence Sx + 2 t|s« = |^[^ 



4S^ 



If now we impose the condition that during the hypothetical 
motion, as during the actual motion, the equation of the con- 
servation of energy holds good, that is, that 

then ST + Sr = 0, 

• SZ = ST-Sr = 2ST = S(2T), 
and (1) becomes 
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and, finally, S C'2Tdt = 0. (2) 

The equation -^ = I 2Tdt defines the action^ A^ of the forces, 

and the fact stated in (2) is usually called the principle of least 
action. As a matter of fact, (2) shows merely that the action is 
" stationary." 



45. In establishing Hamilton's principle we have supposed 
the course followed by the system to be varied, subject merely 
to the limitation that the time of transit from initial to final 
configuration should be unaltered ; consequently, as the total 
energy is not conserved, the varied course is not a natural 
course. That is, to compel the system to follow such a course 
we should have to introduce additional forces that would 
do work. 

In establishing the principle of least action, however, we 
have supposed the course followed by the system to be varied, 
subject to the limitation that the total energy should be unal- 
tered ; consequently the varied course is a natural course. That 
is, to compel the system to follow such a course we need 
introduce merely suitable constraints that would do no work. 

We have deduced both principles from the equations of 
motion. Conversely, the equations of motion can be deduced 
from either of them. Each of them is therefore a necessary 
and sufiicient condition for the equations of motion. 
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46. The action, Ay of a conservative set of forces acting on a 

moving system has been defined as the time integral of twice 

the kinetic energy. ^ 

A= f '2Tdt. (1) 

But A= f ''^m (d^-i-f-h i') dt= C \m (^\ dt= C \mv^dt 



Therefore 



= C^mvds, (2) 



so that the action might just as well have been defined as the 
sum of terms any one of which is the line integral of the momen- 
tum of a particle taken along the actual path of the particle. 

There is another interesting expression for the action, which 
does not involve the time even through a velocity. 

Since T+r=A, 2T=2(A-r). 



m 



As T=^'^(i? + f + i^, 

£.(|)'.2(*-r), 



and 



Y.mds' ■ 


= 2(h-r)dt\ 


dt: 
A: 


~N2(A-r)' 

r . 



A = Cy/2(h- r)^md8\ (3) 

47. We have stated without proof that the differential equa- 
tions of motion for any system under conservative forces could 
be deduced from Hamilton's principle or from the principle of 
stationary action. Instead of giving the proof in general, we will 
give it in a concrete case, that of a projectile under gravity. 
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We shall use fixed rectangular coordinates, F horizontal and 
X vertically downward, taking the origin at the starting point 
of the projectile. In this case, obviously, 

T=|(i^ + y^), and V=-w^. 
(a) By Hamilton's principle, we have 

hp'jii!'^f + 2gx']dt^Q, (1) 

or r'i[f + f^2gx']dt = (). 

I \_xSx -|- ^S^ -|- gSx] dt = 0. 

j\xj^Sx-\-if—Bi/+gBx\dt = 0. 
Integrating by parts, 

but since Sx and Sy are wholly arbitrary, this is impossible unless 
the coefficients of Sx and Sy in the integrand are both zero. 

Hence ^-^ = 0, 

dt' "• 
(6) By the principle of stationary action, we have 

whence C'C^ + f)i^dr-^r'^h(f + f)dr = (i. (1) 
Jr. «»• Jr. ar 
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fit 
But o" (^ "*" ^) "~ ^^^^ ~ ^ 

Hence S ^^ — = gix^ 

and S(i« + ^)=8f±£+^8a;. 

^_l_ ^ d d d 

sm<« ^-^di~di^-Ttdt^' (^- App. B, § 6, (11) 

and (1) becomes 



or 






Integrating by parts, we get 
^xSx + ifSi/y-JY^ 

But 52; and Sy are zero at the beginning and at the end of 
the actual and the hypothetical paths, so that 

X^i dt 

i(x-g)hx + yhy']—dr = 0, 
11 

and as hx and hy are wholly arbitrary, this necessitates that 

y = o, 

as in Art. 47, (a). 
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48. Although we have proyed merely that in a system under 
conservative forces the action satisfies the necessary condition 
of a minimum, namely, SA = 0, it may be shown by an elabo- 
rate investigatioD that in most cases it actually is a minimum, 
and that the name ** least action,'^ usually associated with the 
principle, is justified. 

A very pretty corollary of the principle comes from its applica- 
tion to a system mo\'ing under no forces or under constraining 
forces that do no work. In either of these cases the potential 
energy V is zero, and consequently the kinetic energy is constant, 
thatis,T=t , , 

A= r'2Tdt=2h r'dt = 2h(t^-^Q 

and the action is proportional to the time of transit. Hence 
the actual motion is along the course which occupies the least 
possible time. 

For instance, if a particle is moving under no forces, the 

energy -^ v* is constant, the velocity of the particle is constant, 

and as it moves from start to finish in the least possible time, 
it must move from start to finish by the shortest posi^ible path, 
that is, by a straight line. 

If instead of moving freely the particle is constramed to 
move on a given surface, the same argument proves that it 
must trace a geodetic on the surface. 

49. Varying Action.* The action, A, between two configura- 
tions of a system under conservative forces is theoretically 
expressible in terms of the initial and final coordinates and the 
total energy (v. Art. 46, (3)), and, when so expressed, Hamilton 
called it the characteristic function. 

In like manner, if 

S= f*'Ldt= r\T-V^dt=^ f\f^- 2V}dt, 

* For a detailed account of Hamilton's method, see Kouth, Advanced Rigid 
Dynamics, chap, x, or Webster, Dynamics, § 41. 
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S may be expressed in terms of the initial and final coordinates 
of the system, the total energy, and the time of transit. When 
so expressed, Hamilton called it the principal functioTL By 
considering the variation produced in either of these functions 
by varying the final configuration of the system, Hamilton 
show^ed that from either of these functions the integrals of 
the differential equations of motion for the moving system 
could be obtained, and he discovered a partial differential 
equation of the first order that S would satisfy and one that A 
would satisfy, and so reduced the problem of solving the equa- 
tions of motion for any conservative system to the solution of 
a partial differential equation of the first order. 

It must be confessed, however, that in most cases the advan- 
tage thus gained is theoretical rather than practical, as the 
solving of the equation for S or for A is apt to be at least as 
difficult as the direct solving of the equations of motion. 



CHAPTER V 

APPLICATION TO PHYSICS 

50. Concealed Bodies. In many problems in mechanics the 
configuration of the system is completely known, so that a set 
of coordinates can be chosen that will fix the configuration 
completely at any time; and the forces are given, go that if 
the system is conservative the potential energy can be found. 
In that case the Lagrangian function L or the Hamiltonian 
function H can be formed, and then the problem of the motion 
of the system can be solved completely by forming and solving 
the equations of motion. 

In most problems in physics, however, and in some problems 
in mechanics the state of things is altogether different. It may 
be impossible to know the configuration or the forces in their 
entirety, so that the choice of a complete set of coordinates or 
the accurate forming of the potential function is beyond 'our 
powers, while it is possible to observe, to measure, and partly to 
control the phenomena exhibited by the moving system which 
we are studying. If from results which have been, observed or 
have been deduced from experiment we are able to set up in- 
directly the Lagrangian function, or the Hamiltonian function, 
or the modified Lagrangian function, we can then form our dif- 
ferential equations and use them with confidence and profit. 

51. Take, for instance, the motion considered in the problem 
of the two equal particles and the table with a hole in it (v. 
Art. 8, (<i)), and suppose the investigator placed beneath the 
table and provided with the tools of his trade, but unable to 
see what is going on above the surface of the table. He sees 
the hanging particle and is able to determine its mass, its 

08 
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velocity, and its acceleration ; to fix its position of equilibrium ; 
to measure its motion under various conditions ; to apply addi- 
tional forces, finite or impulsive, and to note their effects. 

His system has apparently one degree of freedom. It, pos- 
sesses an apparent mass m and is certainly acted on by gravity 
with a downward force mg. He determines its position of 
equilibrium, and taking as his coordinate its distance x below 
that position, he painfully and laboriously finds that i, the accel- 
eration of the hanging mass, is equal to 






^• 



dx dx 



He is now ready to call into play his knowledge of mechanics. 
If T is the kinetic energy and L is the Lagrangian function for 
the system which he sees, 

and L=T-V=^x'-V. 

dt\dx/ dx 

Therefore g = |[_^_lj^. 

The motion, then, can be accounted for as due to the 
downward force of gravity combined with a second vertical 

force having the potential energy -^ — -5 + a; , that is, a 

force vertically upward having the intensity -^ -^ "*" ■'^ ' 

and of course this force must be the pull of the string and may 
be due to the action of some concealed set of springs. 
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On the other hand, the moving system may contain some 
concealed body or bodies in motion, and that this is the fact is 
strongly suggested if a downward impulsive force P is applied 
to the hanging particle when at rest in its position of equilib- 
rium. For such a force is found to impart an instantaneous 

P 

velocity x = - — » just half what we should get if the hanging 

particle were the only body in the system, and just what we 
should have if there were a second body of mass m above the 
table, connected with the hanging particle by a stretched string 
of fixed length. 

Obviously this concealed body is igrwrahle^ since we have 
already found a function from which we can obtain the differ- 
ential equation for x^ namely, the function we have just called 
the Lagrangian function X. It is 






and contains the single coordinate x. 

But if we have ignored a concealed moving body form- 
ing part of our system, this expression is not the Lagrangian 
function Z, but is equal or proportional to <E>, the Lagrangian 
function modified for the coordinate or coordinates of the con- 
cealed body ; and in that case some or all of the terms which 
we have regarded as representing the potential energy of the 
system may be due to the kinetic energy of the concealed body 
(v. Art. 41). 

Of course the complete system may have two or more degrees 
of freedom. Let us see what we can do with two. Take x 
and a second coordinate and remember that as ^ is ignor- 
able it must be a ct/clic coordinate and must not enter into 
the potential energy. 

Let us now form the Lagrangian function and modify it for 0. 

We have T — A^ -h Bxd -|- (7^^ where A^ B^ and C are functions 
of X. 
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A_ P« — ^i 

20 ' 

To modify for we must subtract ^p^ We get 

2C^ ^ 4C 2C 



'V-^k 






Suppose that no external force acts on the concealed particle, 
so that the potential energy of the system is — mgx. 
Then, if the Lagrangian function modified for ^ is 4>, 

Since on our ignoration hypothesis 6 does not enter the 
potential energy, the momentum Pq = K^ a constant, and 



* = M-T— i^-hTT^^iJ-TT,"'- '^9^' 



2C 4(7 



But we know that ^ is equal or proportional to 

m ^ mg a' mgx 

the function which on our hypothesis of no concealed bodies 
we called L. 

We see that if -5 = 0, if ^ = tw, and if -j— - = .^ ^ — — , 

4 C/ 2 (a — a;)^ 

whence C = -t — ^^ — j-^ » 

<E) = 2Z. 

Then we have T =^ mi? + )-^-^~^ 6", 

2 mgar 
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where K is the momentum corresponding to the coordinate 
and may be any constant. 

If we take for K^ the value m^ga\ 

is obviously the kinetic energy of a mass m moving in a plane 
and having a — x and as polar coordinates. The Lagrangian 
function is equal to 

and the x Lagrangian equation is 

2 mx = — m (a — x')d^ + mg. (1) 

If the angular velocity of the concealed mass when the hang- 
ing particle is at rest in its position of equilibrium is 6^, since 

in that case x=0 and ^ = 0, equation (1) gives us 6^ = -^- • The 

observed motion of the hanging particle is then accounted for 
completely by the hypothesis that it is attached by a string 
to an equal particle revolving on the table and describing a 
circle of radius a about the hole in the table with angular 

velocity -vj- when the hanging particle is at rest in its posi- 
tion of equilibrium. We see that on this hypothesis the term 
o o7 \2"*"^ ' which on the hypothesis that the system 

contained only the hanging particle was an unforeseen part of 
the potential energy, is due to the kinetic energy of the con- 
cealed moving body. 

It may seem that giving K a different value might lead to 
a different hypothesis as to the motion of the concealed body 
that would account for the motion of the hanging particle. 
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Such, however, is not the case. We have 

Let <f, = -^=0. 

Then r = ^ [ir* + ir* + (a - a:)» i^"]. 

~ [ir* + (a — ^y <f>^'] is, as above, the kinetic energy of a body 

of mass w, with polar coordinates a — x and <^ ; and the concealed 
motion is precisely as before. In using the form (2) we have 
merely used as our second generalized coordinate ^, a perfectly 
suitable parameter but one less simple than the polar angle. 

52. Problems in Physics. In physical problems there may be 
present electrical and magnetic phenomena and concealed 
molecular motions, as well as the visible motions of the mate- 
rial parts of the system. In such cases, to fix the configuration 
of the system even so far as it is capable of being directly 
observed we must employ not only geometrical coordinates 
required to fix the positions of its material constituents, but also 
parameters that will fix its electrical or magnetic state ; and as 
we are rarely sure of the absence of concealed molecular 
motions, we must often allow for the probability that the func- 
tion we are trying to form by the aid of observation and 
experiment and on which we are to base our Lagrangian equa- 
tions of motion may be the Lagrangian function modified 
for the ignored coordinates corresponding to the concealed 
molecular motions. 

53. Suppose, for instance, that we have two similar, parallel, 
straight, conducting wires, through which electric currents due 
to applied electromotive forces are flowing. It is found ex- 
perimentally that the wires attract each other if the currents 
have the same direction, and repel each other if they have 
opposite directions, and that reversing the currents without 
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altering the strength of the applied electromotive forces does 
not affect the observed attraction or repulsion. It is known 
that an electromotive force drives a current against the resist- 
ance of the conductor, and that the intensity of the current is 
proportional to the electromotive force. Moreover, it is known 
that an electromotive force does not directly cause any motion 
of the conductor. It is found that, as far as electric currents 
are concerned, the phenomena depend merely on the intensity 
and direction of the currents. 

To fix our configuration we shall take x as the distance 
between the wires and take parameters to fix the intensities of 
the two currents. These parameters might be regarded as 
coordinates or as generalized velocities, but many experiments 
suggest that they are velocities. We shall call them y^ and y^ 
and define y^ as the number of units of electricity that have 
crossed a right section of the first wire since a given epoch. 

As all effects depend upon ij^ and y^ and not on y^ and y^, 
y^ and y^ are cyclic coordinates. 

Let us suppose that there are no concealed motions. Then 

the kinetic energy T is a homogeneous quadratic in x^ y^^ and 

y^. Let 

T^A^^Lyl-VMyJj^^Nyl^Bxy^^ Cxy^, (1) 

where the coefficients are functions of x. 

Since reversing the directions of the currents, that is, revers- 
ing the signs of y^ and y^^ does not change the other phenomena, 
it must not affect T\ therefore B and C are zero. 

If the wires are not allowed to move, i = and T reduces 
to Lyl + My^y^ -}- Ny^^ which is called the electrokinetic energy 
of the system. Since from considerations of symmetry this can- 
not be altered by interchanging the currents, N= L. 

Let us now suppose that the first wire is fastened in position 
and that the only external forces are the electromotive forces 
E^ and E^, producing the currents in the two wires ; the resist- 
ances Ey^ and Ey^ of the wires, equal, respectively, to U^ and 
U^ when the currents y^ and y^ are steady; and an ordinary 
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mechanical force F^ tending to separate the wires. We now 
have ^rp 

and for our x Lagrangian equation 

c\ i" . ckdA .« dA .« dL .« dM . . dL ,« _ ,^. 
2^^ + 2-^--a:»--y/--^,,^,,--y« = i^. (2) 

Let us study this equation. First suppose the electromotive 
forces and the impressed force F all zero, so that y^^y^ = 0, 
and F=0. Equation (2) reduces to 

2Ax-\-^d? = Q, 
dx 

1 dA .« 

or x=- ir, 

. 2Adx 

If, then, a transverse velocity were impressed on the second 

wire, the wire would have an acceleration unless --- = 0. But 

dx 

both wires, on our hypothesis, being inert, they can neither 
attract nor repel each other, therefore -4 is a constant; and as 
T=Ax^ when y^ = y^ = 0^ A is a positive constant. Therefore 

Let us now suppose that ^2 = 0, and F=0. Equation (3) 

becomes 2 Ax — -- y^* = 0, 

dx 

1 dL .„ 
^ = 21^^- 

and the second wire is attracted or repelled by the first, 
even when no current is flowing through the second wire. 

This is contrary to observation. Therefore -r- = 0, and L is 

dx 

a constant; and as T=Lyl when x and y^ are zero, X is a 
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positive constant. Equation (3) now becomes 

and if 2^=0. -i^^^jj,^. 

If y. and y^ are of the same sign, x and -r— have the same 
* * ax 

sign. But according to observation the wires attract if the 
currents flow in the same direction, therefore x is negative and 

--- IS negative. 
dx ^ 

If 2; = 0, we have a single wire carrying a current of inten- 
sity if^ -h ^j. The electrokinetic energy Ly\ + My^^ + Lyl be- 
comes -^(^1 + ^2)^ *^^ -^y? + 2 X^j^j + X^a , so that the value of 
M when a; = -is 2 Z. Hence, as i(f is a decreasing function, M 
is less than 2 L always. As -If is easily seen to be zero when 
X is infinite, it must be positive for all values of x. 

We have, then, 

T^A^^ Lyl + My^y^ + Lyl (4) 

where A and L are positive constants and ^ is a positive decreas- 
ing function of x always less than 2 Z. Z is called the coefficient 
of self-induction of either wire per unit of length, and M the co^- 
dent of mutual induction of the pair of wires per unit of length. 
Our Lagrangian equations are 

'^^i-^vA-'F, (5) 

|[2Xy, + JfyJ = JF,-i?^,, (6) 

jPLy^+My;\ = E^-Ry^. (7) 

S4. Induced Currents, (a) Suppose that no current is flow- 
ing in either of the wires considered in the last section, and 
that the first wire is suddenly connected with a battery fur- 
nishing an electromotive force jE'^, and that thereby a current 
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if^ is impulsively established. Then, by Thomson's Theorem 
(v. Art. 30), such impulsive velocities must be set up in the 
system as to make the energy have the least possible value 
consistent with the velocity caused by the applied impulse. 
If the current y^ set up in the first wire has the intensity i^ 

and making T a minimum, we have 

2^i = 0, 
Mi+2Ly^ = 0; 

whence f = 0, and the second wire will have no initial velocity, 

Mi 

and ^2 = — -r— > and a current will be set up impulsively in the 

second wire, of intensity proportional to the intensity of the cur- 
rent in the first wire. Since, as we have seen, M and L are both 
positive, y^ is negative, and this so-called induced current will be 
opposite in direction to the impressed current y^. This impulsively 
induced current is soon destroyed by the resistance of the wire. 
(6) Suppose that a steady current y^^ caused by the electro- 
motive force -&j, is flowing in the first wire while the second 
wire is inert, and that E^^ and consequently y^^ is impulsively 
destroyed by suddenly disconnecting the battery. This amounts 
to impulsively applying to the first wire the additional electro- 
motive force —E^. If the system were initially inert, this, as 

we have just seen in (a), would immediately set up the induced 

Mi , 
current y^ = ~ — in the second wire, and we should have y^ = — f, 

Mi . 

y^ = -T — , as the immediate result of our impulsive action. Combine 

this with the initial motion in our actual problem, y^ = i, y^ = 0, 

Mi 
and we get for our actual result y^ = 0, y^ = -— : (v. Art. 34). 

So that if our first wire is suddenly disconnected from the 
battery, an induced current whose direction is the same as that 
of the original current is set up in the second wire. It is, how- 
ever, soon destroyed by the resistance of the wire. 
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((?) Suppose that we have a current y^ in our fixed wke, caused 
by a battery of electromotive force E^^ and no current in our 
second wire, and that the second wire is made to move away 
from the first. Equations (6) and (7) of the preceding section 
give us 



dt 



-(2Ly,-My;)i^, 



(4 z« - Jir") ^ = 2 i (je; - Ey;)- M^E^ - By;) 

-(^iLy^-My;ix^£. 

When we are starting to move the Becond wne, 

-^", = 0, ^, = 0, E^-Ry^^Q, 
and we have 

(4 L^ _ Jf ») ^ = My^x ^, 

Civ Ol/X 

As we have seen, Z, Jtf", 4 L^ — -M^, are positive and —— is negar 

ax 

tive. Hence the current y^ will decrease in intensity, and a 
current y^ having the same direction as y^ will be induced in 
the moving wire. 

The phenomena of induced currents which we have just 
inferred from our Lagrangian equations are entirely confirmed 
by observation and experiment. 
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SYLLABUS. DYNAMICS OF A RIGID BODY 

1. D'Alembert's Principle. In a moving system of particles the 
resultant of all the forces external and internal that act on any 
particle is called the effective force on that particle. Its rectangular 
components are rrdi, my, and mz. 

The science of rigid dynamics is based on lyAlewher^s prin- 
ciple : In any moving system the actual forces impressed and 
internal, and the effective forces reversed in direction, form a set 
of forces in equilibrium, and if the system is a single rigid body, 
the internal forces are a set separately in equilibrium and may 
be disregarded. 

It follows from this principle that in any moving system the 
actual forces and the effective forces are mechanically equivalent. 
Hence, 

(a) Smd; = 2X, 

{b) Sm lyx - iri/] = 2 iyX - xY], 

(c) 2m [aJScc + yhj + «8«] = 2 [J^&c + FSy + Zhz], 

These equations may be put into words as follows : 

(a) The sum of those components which have a given direction 
is the same for the effective forces and for the actual forces. 

(h) The sum of the moments about any fixed line is the same for 
the effective forces and for the actual forces. 

(c) In any displacement of the system, actual or hypothetical, the 
work done by the effective forces is equal to the work done by the 
actual forces. 

Equations (a) and (h) are called differential equations of motion 
for the system. 

2, p^ = ^mv^ = Imx and is the linear mxymentum, of the system in 
the X direction ; h^ = 2m [yv^ — xv^"] = 2m [jjx — xj/"] and is the 
moment of momentum, about the axis of Zm 

109 
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Equations (a) and (b) of § 1 may be written, respectively, 

f = 5Z, and f = S[y^-.y]. 

and § 1, (a), and § 1, (6), may now be stated as follows : 

(a) In a moving system the rate of change of the linear momentum 
in any given direction is equal to the sum of those components of 
the actual forces which have the direction in question. 

(b) In a moving system the rate of change of the moment of 
momentum about any line fixed in space is equal to the sum of the 
moments of the actual forces about that line. 

3. Center of Gravity. x=^—-—' 

M 

Hence v^ = Smi = ilf -- ; 

at 

or, the linear momentum in the X direction is what it would be if 
the whole system were concentrated at its center of gravity. 

or, the moment of momentum, about the axis of Z is what it would be if 
the whole mass were concentrated at the center of gravity plus what 
it would be if the center of gravity were at rest at the origin and the 
actual motion were what the relative motion about the moving center 
of gravity really is. 

4. The motion of the center of gravity of a moving system is the 
same as if all the mass were concentrated there and all the actual 
forces, unchanged in direction and magnitude, were applied there. 

The motion about the center of gravity is the same as if the center 
of gravity were fixed in space and the actual forces were unchanged 
in magnitude, direction, and point of application. 

5. If the system is a rigid body containing a fixed axis, 

where Mk^ = M(h^ + k^) and is the moment of inertia about the axis, 
and where w is the angular velocity of the body. 

Equation (b) of § 1 becomes Mk*^ — = iV, where N is the sum 

of the moments of the impressed forces about the fixed axis. 
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6. If the system is a rigid body containing a fixed point and o),, 
a>y, Wg, are its angular velocities about three axes fixed in space and 
passing through the fixed point, 

where C is the moment of inertia about the axis of Z, and D and E 

are the products of inertia about the axes of X and y, respectively ; 

that is, 

C = Sm (a:^ + y^, D = ^myz, E = ^Sanzx, 

Equation (h) of § 1 becomes 

dm. 



dh dta^ dm^ 



^__ 



(^ — B) <Og.iOj, + E<aj,<Dg 



dt dt dt dt 

where N is the sum of the moments of the impressed forces about 
the axis of Z. 

7. Euler's Equations. If the system is a rigid body containing 
a fixed point and w^, w^, Wg, are its angular velocities about the 
principal axes of inertia through that point (a set of axes fixed in 
the body and moving with it), equation (h) of § 1 becomes 

8. Euler's Angles. Euler's angles 6, ij/, <f>y are coordinates of a mov- 
ing system of rectangular axes A, B, C, referred to a fixed system 
-Y, F, Z, having the same 
origin 0. 9 is the colati- 
tude and ij/ the longitude 
of the moving axis of C 
in the fixed system (re- 
garded as a spherical sys- 
tem with the fixed axis of 
Z as the polar axis), and 
<!> is the angle made by 
the moving C-4-plane with 
the plane through the fixed 
axis of Z and the moving 
axis of C 
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\\ e hare ^ ^ ^ ^^^ ^ ~~ ^ ^^ ^ <^<^ ^7 

«^ = ^ cos ^ + ^ sin $ sin ^^ 
•i^ = ^C08tf + ^; 
and M^ = ~ ^ sin ^ + ^ sin tf cos ^, 

M, = tf COS ^ + ^ sin tf sin ^^ 
•, = ^ cos * -h ^. 

9. 4^ wi**, the kinetic energy of a particle y becomes 
7, [-r^ + y* -h 2^ in rectangular coordinates, 

MA 

— [r^ + r*^*] in polar coordinates, 

ad 

-^ [r^ + 1^1^^ + sin-^*)] in spherical coordinates. 

10. 2 "o" '^J ^^^ kinetic energy of a moving system y becomes 

V-^(i* + y^ + i^in rectangular coordinates ; 
^ Ar^w* if a rigid body contains a fixed axis ; 

V (~//) '^Vif) I "^ 'o^^ ^ *^® body is free and the motion 
is two-dimensional ; 

|[.4a>/+5«/+C«.«-22>«,«.-2£ai.«, -2Fai,a)J if the 

body is rigid and contains a fixed point ; 

^[^«>i*+5<u2* + Ca%*] if the body is rigid and the axes are 
the principal axes for the fixed point. 

11. ImpnlsiTe Forces. In a system acted on by impulsive forces, 
the resultant of all the impulsive forces external and internal that 
act on any particle is called the effective impulsive force on that 
jjarticle. Its rectangular components are m (x^ — x^), m (y^ — y^), 
(?fiz^ — z^). D'Alembert's principle holds for impulsive forces. It 
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follows that in any system acted on by imptilsive forces, the actual 
forces and the effective forces are mechanically equivalent. HencOi 

(a) 2m [i;^ - X J = 2X, 

(b) Sm [y^x^ - x^y^ - y^x^ -f x^^"] = 2 [yX - x F], 

(c) 2m l(x^ - x^) &c 4- (y^ - y^) 8^ + (k^ - ;^^) 82] 

= 2[-X:&c4-F8y+Z8«]. 

These equations may be put into words as follows : 

(a) The sum of the components which have a given direction 
is the same for the effective impulsive forces and for the actual 
impulsive forces. 

(b) The sum of the moments about a given line is the same for 
the effective impulsive forces and for the actual impulsive forces. 

(c) In any displacement of the system, actual or hypothetical, 
the sum of the virtual moments of the effective impulsive forces is 
equal to the sum of the virtual moments of the actual impulsive 
forces. 

Equations (a) and (b) are the equations for the initial motion 
under impulsive forces, (a) and (b) may be restated as follows: 

In a system acted on by impulsive forces, the total change in the 
linear momentum in any given direction is equal to the sum of those 
components of the actual impulsive forces which have the direction 
in question. 

In a system acted on by impulsive forces, the total change in the 
moment of momentum about any fixed line is equal to the sum of 
the moments of the actual impulsive forces about that line. 

Section 4 holds unaltered for impulsive forces. 
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THE CALCULUS OF VARIATIONS 

1. The calculus of variations owed its origin to the attempt to 
solve a very interesting class of problems in maxima and minima 
in which it is required to find the form of a function such that the 
definite integral of an expression involving that function and its 
derivatives shall be a maximum or a minimum. 

Take a simple case: If y=^f(x\ let it be required to deter- 
mine the form of the function /, so that I <^ aj, y, -^ L/a; shall be 
a maximum or a minimum. *^'o 

Let f{x) and F(x) be two 
possible forms of the function. 
Consider their graphs y =f(x) 
and y = F(x). 

Urj(a')is F(x)^f(x),rf(x) 
can be regarded as the incre- 
ment given to y by changing 
the form of the function from 
f(x) to F(x)y the value of the 
independent variable x being 
held fast. 

This increment of y, rj (x), is called the variation of y and is written 

Sy ; it is a function of x, and usually a wholly arbitrary function of x. 

dv 
The corresponding increment in y\ where y' = -f^ > can be shown to be 

dv 
rj'(x)f and is the variation of y\ and is written Sy', ot B-f-' Obviously, 




'y'=fjy- 



(1) 



If an infinitesimal increment hy is given to y^ it is proved in the 

differential calculus that -7- <f>(y)^y differs by an infinitesimal of 

higher order than Sy from the increment produced in <f>(y)* This 
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approximate increment is called the variation of <^ (2^), so that 

8*(y) = |;«^(2/)8y- (2) 

Similarly, 8,^(y, y')=^8y + ^,8y', (3) 

or since x is not varied, 

8<^(^,y,y')=^8y + |^,8y'. (4) 

As d<^ (y) = — <^ (y) dyy 

and d^ {y, y') = -^dy + '^, dy\ 

we can calculate variations by the familiar formulas and processes 
used in calculating differentials. 

2. Let a be an independent parameter. Then y =if(x)-\- (Xrj(x), 
or y =f(x) + aSy, is any one of a family of curves including 
y =f(x) (corresponding to a = 0) and y =f(x) + rj (x) (corresponding 
to a: = 1). If Xq and x^ are fixed values, and if 

I <^ (ic, 3/ -f a^y, y' + aSy') dx, 

1(a) is a function of the parameter a only. A necessary condition 
that /(a), a function of a single variable er, should be a maximum 
or a minimum when a = is /'(a) = when a = 0. 



'«-X"K''^v'^']^ 



when a = 0. 

That is, ^'(0)= f '3<^rfa;. (v. § 1, (4)) 

l <l>(x, y, 2/')^^ should be a 



maximum or a minimum 



when y =f(x) is I 8<^(a;, y, y')dx = 0. 
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j if^dx is taken as the definition of the variation of / <l>(ixy and 
our necessary condition is usually written 

How this condition helps toward the determination of the form 
of /(jr) can be seen from an example. 

3. Let it be required to find the form of the shortest curve 
y =/(a;) joining two given points {x^, y^ and (x^, y^). Here, since 



-0 

and / is to be made a minimum. 



£& = Vl -f y^'dx, 
I = f'Vl -H y'^dx, 



3/= C sVl + Z'-eia; 

-L 










. 1 dx 

Vl + 2/'* 

— cfe 



Integrating by parts, this last reduces to 



since, as the ends of the path are given, 8y = when x^x^ and 



when X = x^. Then 8/ = if 



• ^ 82/rfx = 0; 






c?a; Vl -f- y'^ 
but since our Sy (that is, i; (a)), is a function which is wholly arbitrary, 
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d v' 
the other factor, , must be equal to zero if the integral 

is to vanish. '^ ^-r y 

7/' 

This gives us , ^ = C 

vrT7 

and the required curve is a straight line. 

4. In our more general problem it may be shown in like manner that 






leads to a differential equation between y and x and so determines y 
as a function of x. 

Of course 8/ = is not a sufficient condition for the existence of 
either a maximum or a minimum and does not enable us to discrim- 
inate between maxima and minima, but like the necessary condition 

-~ = for a maximum or minimum value in a function of a single 

variable, it often is enough to lead us to the solution of the problem. 

5. Let us now generalize a little. Let a;, 3^, ;sj, • • ., be functions of 

dx d\i d% 

an independent variable r, and let sc' = — > y* = -j-y «' = —-,.... 

Suppose we have a function <^(r, a, 2^, «, • . ., x\ y\ z\ • • •). 

By changing the forms of the fimctions but holding r fast, let 
^j 3^> ^? • • •> ^ given the increments ^(r), 17 (r), ^(r), • • .. 

X then becomes x -f- f (r), y becomes y + 17 (r), z becomes « + ^ (r), • • • ; 
sc' becomes aj' + f '(r), 3/' becomes 3/' -f- 17' (r), «' becomes »' 4- ^W? • • •• 

i(r), $'(r), are the variations of x and aj' and are written &c 

and &c'. Obviously, d 

&k' = -- &c. 
ar 

The increment produced in <^ when infinitesimal increments &r, 
81/, • . ., &c', 8y', . . ., are given to the dependent variables and to their 
derivatives with respect to r is known to differ from 

by terms of higher order than the variations involved. This approxi- 
mate increment is called the variation of if> and is written 8<^. It is 
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found in any case precisely as c?^, the complete differential of ^, 
is found. 

That I 8^(r, «, y, • • •, x', y', • • • ) c?r = is a necessary condition 

that i ^(r, aj, y, • • •, x', y', ' > •) dr should be a maximum or a 

minimum can be established by the reasoning used in the case of 

^(cc, y, y^)dx. The integral I 8<^c?r is called the variation 

I ^c?r, so that 3 I ^c?r = I B<l>dr. 
6. It should be noted that our important formulas 

r.dX __^CS 

dr dr ' 
and 3 I <l}dr= I h<f>dry 



hold only when r is the independent variable which is held fast 
when the forms of the functions are varied; that is, when 8r is 
supposed to be zero. , 

If X and y are functions of r and we need 8-^> we get it in- 
directly thus: ^ f 

dx cc' 

"^ = ^^ = ""^^^ "^^'^^ -'xid'r^^~'d^Try 

so that ^3^ = 3-^^-T^T-^- (1) 

aic dx dx dx ^ 

If we need 8 I ^6?ic, we must change our variable of integration 

8 I <l>dx=:S j 4t—dr. 
^Ji,dx=jh{^^^dr. (2) 
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